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1 Preface 

The dedicated study of 6-flavoured hadrons has developed into one of the most active and 
most promising areas of experimental high-energy physics. The detailed investigation of b 
decays at the 5-meson factories and at hadron colliders will probe the flavour physics of 
quarks with unprecedented precision. To fully exploit this rich source of data a systematic 
theoretical approach is necessary. The required field theoretical tools are the subject of 
these lectures. 

We shall discuss the construction of effective weak Hamiltonians, introducing the op- 
erator product expansion (OPE) and the renormalization group (RG) and presenting the 
effective Hamiltonians for nonleptonic AB = 1 and AB = 2 transitions as examples. The 
subsequent chapter is devoted to heavy-quark effective theory (HQET). It explains the 
basic formalism as well as the application to heavy-light and heavy-heavy currents, dis- 
cussing the 5-meson decay constant /# and the exclusive semileptonic decay B — > D*lu, 
respectively. Inclusive b decays and the heavy-quark expansion (HQE) are treated next, in 
particular the general formalism, the issue of quark-hadron duality, the theory of 6-hadron 
lifetimes and of inclusive semileptonic decays. Finally, we discuss QCD factorization for 
exclusive hadronic B decays, focussing on B — > Dn and B — > nn. We conclude with a 
short summary. 

The effective-Hamiltonian framework is the oldest and most general of the methods 
we shall discuss. It dates back, more or less, to the beginnings of the standard model 
itself. HQET and HQE are later developments that have been established in the second 
half of the eighties and at the beginning of the nineties. The last topic, QCD factorization 
for exclusive B decays is the most recent. It is the least well established among these 
methods and it continues to be studied and developed in further detail. 

We would like to mention a very short selection of literature, which we hope will be 
helpful to obtain further information on the various subjects related to the contents of 
these lectures. Very useful resources are the BaBar Physics Book (Harrison & Quinn 
1998) and the Fermilab B Physics Report (Anikeev et al. 2002). They collect nice 
reviews on both theoretical and experimental topics in B physics. A textbook more 
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specifically directed towards theoretical heavy quark physics is the work of Manohar & 
Wise (2000). Review articles on particular subjects are (Buchalla, Buras & Lautenbacher 
1996) on effective weak Hamiltonians, (Neubert 1994) on HQET and (Bigi et al. 1994; 
Bigi, Shifman and Uraltsev 1997) on HQE. 



2 Introduction and overview 



2.1 Motivation 



In the following chapters we will study the theoretical tools to compute weak decay 
properties of heavy hadrons. To put the formalism into perspective, we start by recalling 
the main motivation for this subject. 

The central goal is the investigation of flavour physics, the most complicated sector 
in our understanding of fundamental interactions. A good example is give by particle- 
antiparticle mixing, as first studied with neutral kaons. The strong interaction eigenstates 
K° (sd) and K° (ds) can be transformed into each other through second order weak 
interactions, which leads to a tiny off-diagonal entry M\2 in the mass matrix (Fig. [I]). The 
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Figure 1. K - K mixing. 

physical eigenstates are K L S = (K° ± K°)/^/2. Their mass difference Am^ = — ms 
is given by 2|M 12 | and reads 

— « ^ft^ WM 2 ml = 7 • 10- 15 (1) 

171k 07T z 

where the number on the right-hand side is the experimental value. The theoretical 
expression is derived neglecting the third generation of quarks and CP violation, which 
is a valid approximation for Aitlk- The factors fj^Bx [Bk ~ 1 for the purpose of a first 
estimate) account for the binding of the quarks into mesons. A crucial feature of ([[]) is 
the Glashow-Iliopoulos-Maiani (GIM) cancellation mechanism: The orthogonality of the 
quark mixing matrix V^- ii = u, c; j = d, s) leads to a cancellation among the various 
contributions with intermediate up and charm quarks (symbolically the amplitude has 
the form (mm) — (uc) — (cu) + (cc)). For m u = m c this cancellation would be complete, 
giving Arrix = 0. Still, even for m„ 7^ m c , the contributions from virtual momenta 
k ^> m c , m u cancel since both m u and m c are negligible in this case. What is left is a 
characteristic effect proportional to m^, the up-quark contribution being subleading for 
m u Aqcd <S ™>c- This circumstance allowed Gaillard and Lee in 1974 to correctly 
estimate the charm-quark mass m c « 1.5 GeV, before charm was eventually discovered 
in the Fall of the same year. 
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In a similar way the discovery of Bj - Bd mixing by the ARGUS collaboration (Albrecht 
et al. 1987) proved to be another milestone in flavour physics. In full analogy to K — K 
mixing we have 
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(2) 



where now the top-quark contribution is completely dominant. The unexpectedly large 
value observed by ARGUS provided the first evidence that the top-quark mass (m tiPO ; e = 
176 GeV) was comparable to the weak scale and in any case much heavier than anticipated 
at the time. 

These examples show very nicely the "flavour" of flavour physics: Precision observables 
are sensitive probes of high-energy scales and yield crucial insights into the fundamental 
structure of weak interactions. At the same time we see that hadronic effects manifest in 
quantities such as f B , B B , and strong interactions of the participating quarks in general, 
play an important role. Their understanding is necessary to reveal the underlying flavour 
dynamics and is the main subject of heavy quark theory. 

B - B mixing, CP violation in B decays and other loop- induced reactions of 6-flavoured 
hadrons are of great interest and continue to be pursued by numerous experiments. A 
central target is the Cabibbo-Kobayashi-Maskawa (CKM) matrix 
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which parametrizes charged- current weak interactions in the standard model. The second 
equality in (|3|) introduces the approximate form of the Wolfenstein parametrization, where 
the four independent CKM quantities are A, A, g and r]. The unitarity triangle, as defined 
in terms of g and r\ is shown in Fig. 0, indicating the CP violating angles a, (3 and 7. 




Figure 2. Definition of the unitarity triangle. 

The status of the unitarity triangle (in terms of g = g(l — A 2 /2) and fj = r](l — A 2 /2)) 
is displayed in Fig. ^, with input from CP violation in the kaon sector (ek) and from B 
decays (\V ub /V cb \, Am Bd , Am Bs and sin 2(3). 
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Figure 3. Global fit of the unitarity triangle (darker shaded region), without sin 2/3 from 
CP violation in B — > J/^/Ks- The constraint from the world average sin 2/3 (hatched 
areas) is overlaid (Hocker et al. 2001). 

2.2 B decays — overview 

A vast number of different 5-decay observables is available to further test the standard 
scenario of flavour dynamics. One may distinguish the following broad classes. 

• dominant decays 

b — > cud, b — ► ccs, b — > civ 

B -> Dn, B -> VK, B -> D w /z/ 

• rare decays 

6 — ► MMff, 6 — > mts, b — > -u/z/ 

B -> 7T7T, 5 -> TrK, 5 -> 7T/P, 5 -> Zl? 

• rare and radiative (loop induced) decays 
b — > s(d)j, b — > s(d)l + l~, b — > s{d)vv 

B -> K«7, 5 -> p7, . . . 

• = 2 oscillations 

B d - B d , B s - B s mixing 

Other obvious classifications are between inclusive and exclusive processes or between 
hadronic and (semi)leptonic decays. 

A few key properties of 6-hadrons enhance considerably the possibilities in B physics 
both experimentally and theoretically. First, the smallness of V c b = 0.04 leads to a long 
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lifetime of tb ~ 1.6 ps. In addition the 6-quark mass is large compared with the QCD 
scale 

m h > A Q cd ~ 0.3GeV (4) 

The exact value of mj depends on the definition. In particular, the running MS mass 
mb(rhb) = 4.2±0.1GeV, the pole mass m^poie ~ 4.8GeV, whereas the mass of the lightest b- 
hadron is m B = 5.28GeV. The smallness of kQ C D/ m b provides us with a useful expansion 
parameter. Together with the property of asymptotic freedom of QCD and a s (m&) <C 1, 
this opens the possibility of systematic approximations, which are exploited in the various 
applications of heavy quark theory. 

3 Effective weak Hamiltonians 

The task of computing weak decays of hadrons represents a complicated problem in quan- 
tum field theory. Two typical cases, the first-order nonleptonic process B° — > 7r + 7r~, and 
the loop-induced, second-order weak transition B~ — > K~vv are illustrated in Fig. |j. The 
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Figure 4. QCD effects in weak decays. 

dynamics of the decays is determined by a nontrivial interplay of strong and electroweak 
forces, which is characterized by several energy scales of very different magnitude, the 
W mass, the various quark masses and the QCD scale: m t , M w 3> m&, m c 3> Aq C d 
m u , m&, (m s ). While it is usually sufficient to treat electroweak interactions to lowest 
nonvanishing order in perturbation theory, it is necessary to consider all orders in QCD. 
Asymptotic freedom still allows us to compute the effect of strong interactions at short 
distances perturbatively. However, since the participating hadrons are bound states with 
light quarks, confined inside the hadron by long-distance dynamics, it is clear that also 
nonperturbative QCD interactions enter the decay process in an essential way. 

To deal with this situation, we need a method to disentangle long- and short- distance 
contributions to the decay amplitude in a systematic fashion. A basic tool for this purpose 
is provided by the operator product expansion (OPE). 

3.1 Operator product expansion 

We will now discuss the basic concepts of the OPE for B meson decay amplitudes. These 
concepts are of crucial importance for the theory of weak decay processes, not only in 
the case of B mesons, but also for kaons, mesons with charm, light or heavy baryons and 
weakly decaying hadrons in general. Consider, for instance, the basic W^-boson exchange 
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Figure 5. OPE for weak decays. 



process shown on the left-hand side of Fig. [5[ This diagram mediates the decay of a b 
quark and triggers the nonleptonic decay of a B meson such as B° — ► 7r + 7r~. The quark- 
level transition shown is understood to be dressed with QCD interactions of all kinds, 
including the binding of the quarks into the mesons. To simplify this problem, we may 
look for a suitable expansion parameter, as we are used to do in theoretical physics. Here, 
a key feature is provided by the fact that the W mass M\y is very much heavier than 
the other momentum scales p in the problem (m&, h-QCDi m u, m d, We can therefore 
expand the full amplitude A, schematically, as follows 
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which is sketched in Fig. ||. Up to negligible power corrections of 0(p 2 /M$y), the full 
amplitude on the left-hand side is written as the matrix element of a local four-quark 
operator Q, multiplied by a Wilson coefficient C. This expansion in 1/Mw is called a 
(short- distance) operator product expansion because the nonlocal product of two bilinear 
quark-current operators (du) and (ub) that interact via W exchange, is expanded into a 
series of local operators. Physically, the expansion in Fig. [5| means that the exchange 
of the very heavy W boson can be approximated by a point-like four-quark interaction. 
With this picture the formal terminology of the OPE can be expressed in a more intuitive 
language by interpreting the local four-quark operator as a four-quark interaction vertex 
and the Wilson coefficient as the corresponding coupling constant. Together they define 
an effective Hamiltonian TC e ff = C -Q, describing weak interactions of light quarks at low 
energies. Ignoring QCD the OPE reads explicitly (in momentum space) 

Cp 1 % — 

A = -yKdKb p _ M ^ (du) v ^ A (ub) v _ A 



-i-kVZyj3 ■ (Q) + O | 1 - T ) (6) 



with C = 1, Q = (du)v-A(ub)v-A and 



G 

Heff = ^ V ^ V ^( du )v-A(ub) v ^ A (7) 



As we will demonstrate in more detail below after including QCD effects, the most 
important property of the OPE in @ is the factorization of long- and short- distance 
contributions: All effects of QCD interactions above some factorization scale p, (short 
distances) are contained in the Wilson coefficient C. All the low-energy contributions 
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below fi (long distances) are collected into the matrix elements of local operators (Q). 
In this way the short- distance part of the amplitude can be systematically extracted and 
calculated in perturbation theory. The problem to evaluate the matrix elements of local 
operators between hadron states remains. This task requires in general nonperturbative 
techniques, as for example lattice QCD or QCD sum rules, but it is considerably simpler 
than the original problem of the full standard-model amplitude. In some cases also the 
approximate flavour symmetries of QCD (isospin, SU(3)) can help to determine the non- 
perturbative input. This is true in general for hadronic weak decays. A decisive advantage 
of heavy hadrons is the fact that the heavy-quark mass itself is still large in comparison 
to Aqcd- The limit <C 1 can then be exploited, which is achieved, depending 

on the application, by using heavy quark effective theory, heavy quark expansion or QCD 
factorization for exclusive nonleptonic decays. 

The short-distance OPE that we have described, the resulting effective Hamiltonian, 
and the factorization property are fundamental for the theory of B decays. However, the 
concept of factorization of long- and short-distance contributions reaches far beyond these 
applications. In fact, the idea of factorization, in various forms and generalizations, is 
the key to essentially all applications of perturbative QCD, including the important areas 
of deep-inelastic scattering and jet or lepton pair production in hadron-hadron collisions. 
The reason is the same in all cases: Perturbative QCD is a theory of quarks and gluons, 
but those never appear in isolation and are always bound inside hadrons. Nonperturbative 
dynamics is therefore always relevant to some extent in hadronic reactions, even if these 
occur at very high energy or with a large intrinsic mass scale. Thus, before perturbation 
theory can be applied, nonperturbative input has to be isolated in a systematic way, and 
this is achieved by establishing the property of factorization. It turns out that the weak 
effective Hamiltonian for nonleptonic B decays provides a nice example to demonstrate 
the general idea of factorization in simple and explicit terms. 

We would therefore like to discuss the OPE for B decays in more detail, including the 
effects of QCD, and illustrate the calculation of the Wilson coefficients. A diagrammatic 
representation for the OPE is shown in Fig. |6|. The key to calculating the coefficients 



Cj is again the property of factorization. Since factorization implies the separation of 
all long-distance sensitive features of the amplitude into the matrix elements of (Qi), the 
short-distance quantities Ci are, in particular, independent of the external states. This 
means that the Cj are always the same, no matter whether we consider the actual physical 
amplitude where the quarks are bound inside mesons, or any other, unphysical amplitude 
with on-shell or even off-shell external quark lines. Thus, even though we are ultimately 
interested in, e.g., B — > mi amplitudes, for the perturbative evaluation of Cj we are free to 
choose any treatment of the external quarks according to our calculational convenience. 
A convenient choice that we will use below is to take all external quarks massless and 
with the same off-shell momentum p (p 2 ^ 0). 

The computation of the Ci in perturbation theory then proceeds in the following steps: 




Figure 6. Calculation of Wilson coefficients of the OPE. 
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• Compute the amplitude A in the full theory (with W propagator) for arbitrary 
external states. 

• Compute the matrix elements (Qi) with the same treatment of external states. 

• Extract the Cj from A = Ci (Qi). 

We remark that with the off-shell momenta p for the quark lines the amplitude is even 
gauge dependent and clearly unphysical. However, this dependence is identical for A and 
(Qi) and drops out in the coefficients. The actual calculation is most easily performed 
in Feynman gauge. To 0(a s ) there are four relevant diagrams, the one shown in Fig. || 
together with the remaining three possibilities to connect the two quark lines with a gluon. 
Gluon corrections to either of these quark currents need not be considered, they are the 
same on both sides of the OPE and drop out in the Ci. The operators that appear on 
the right-hand side follow from the actual calculations. Without QCD corrections there 
is only one operator of dimension 6 

Qi = (diUi)v-A(ujbj)v-A (8) 

where the colour indices have been made explicit. To 0(a s ) QCD generates another 
operator 

Q 2 = (d],Uj)v-A(Ujbi)y-. A (9) 

which has the same Dirac and flavour structure, but a different colour form. Its origin is 
illustrated in Fig. [7|, where we recall the useful identity for SU(N) Gell-Mann matrices 




Figure 7. QCD correction with colour assignment. 



(diT&uMujTfii) = -^(diudiujbj) + ^(d^-)^) (10) 
It is convenient to employ a different operator basis, defining 

Q± = 9i±9i 

The corresponding coefficients are then given by 

C± = C l ±C 2 (12) 

If we denote by S± the spinor expressions that correspond to the operators Q± (in other 
words: the tree-level matrix elements of Q±), the full amplitude can be written as 

A=(l + 7+ « s In S+ + ( 1 + 7 -« s In ^) S„ (13) 
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Here we have focused on the logarithmic terms and dropped a constant contribution (of 
order oe s , but nonlogarithmic) . Further, p 2 is the virtuality of the quarks and 7± are 
numbers that we will specify later on. We next compute the matrix elements of the 
operators in the effective theory, using the same approximations, and find 

(Q±) = U + l±a s + In q^) ) S± (14) 

The divergence that appears in this case has been regulated in dimensional regularization 
(D = 4 — 2e dimensions). Requiring 

A = C + (Q+> + C_(Q_> (15) 

we obtain 

/If 2 

C ± = l + 7 ±a,lnHr ( 16 ) 

where the divergence has been subtracted in the minimal subtraction scheme. The effec- 
tive Hamiltonian we have been looking for then reads 

n eff = ^=V* d V ub (C + (n)Q + + C-{p)QJ) (17) 



with the coefficients C± determined in ( jig) to C(a s log) in perturbation theory. The 
following points are worth noting: 



• The 1/e (ultraviolet) divergence in the effective theory (|14D reflects the M\y — > oo 
limit. This can be seen from the amplitude in the full theory ([13]), which is finite, but 
develops a logarithmic singularity in this limit. Consequently, the renormalization 
in the effective theory is directly linked to the In dependence of the decay 
amplitude. 

• We observe that although A and (Q±) both depend on the long-distance properties 
of the external states (through p 2 ), this dependence has dropped out in C±. Here we 
see explicitly how factorization is realized. Technically, to 0(oe s log), factorization 
is equivalent to splitting the logarithm of the full amplitude according to 

M 2 M 2 ii 2 

m= = ln^f + ln^ (18) 
— p z fi z —p A 

Ultimately the logarithms stem from loop momentum integrations and the range of 
large momenta, between My/ and the factorization scale \x, is indeed separated into 
the Wilson coefficients. 

• To obtain a decay amplitude from H. e ff in (0)' ^ ne ma t rrx elements (f\Q±\B)([i) 
have to be taken, normalized at a scale /i. An appropriate value for \i is close to 
the 6-quark mass scale in order not to introduce an unnaturally large scale into the 
calculation of (Q). 

• The factorization scale \i is unphysical. It cancels between Wilson coefficient and 
hadronic matrix element, to a given order in a s , to yield a scale independent decay 
amplitude. The mechanism of this cancellation to 0(a s ) is clear from the above 
example (0) - 
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• In the construction of Ti e ff the H^-boson is said to be "integrated out", that is, 
removed from the effective theory as an explicit degree of freedom. Its effect is still 
implicitly contained in the Wilson coefficients. The extraction of these coeffcients 
is often called a "matching calculation" , matching the full to the effective theory by 
"adjusting" the couplings C±. 

• If we go beyond the leading logarithmic approximation 0(a s \og) and include the 
finite corrections of 0(a s ) in (|13"D, (fl4|), an ambiguity arises when renormalizing 
the divergence in (|T4]) (or, equivalently, in the Wilson coefficients C±). This ambi- 
guity consists in what part of the full (non-logarithmic) 0(a s ) term is attributed 
to the matrix elements, and what part to the Wilson coefficients. In other words, 
coefficients and matrix elements become scheme dependent, that is, dependent on 
the renormalization scheme, beyond the leading logarithmic approximation. The 
scheme dependence is unphysical and cancels in the product of coefficients and ma- 
trix elements. Of course, both quantities have to be evaluated in the same scheme to 
obtain a consistent result. The renormalization scheme is determined in particular 
by the subtraction constants (minimal or non- minimal subtraction of 1/e poles), and 
also by the definition of 75 used in D 7^ 4 dimensions in the context of dimensional 
regularization. 

• Finally, the effective Hamiltonian ( |TTD can be considered as a modern version of the 
old Fermi theory for weak interactions. It is a systematic low-energy approximation 
to the standard model for 6-hadron decays and provides the basis for any further 
analysis. 



3.2 Renormalization group 

Let us have a closer look at the Wilson coefficents, which read explicitly 



where we have now specified the exact form of the 0(a s \og) correction. Numerically 
the factor a s {mb)^± 1 / '(8tt) is about +3.5% (—7%), a reasonable size for a perturbative 
correction (we used a s (fi = 4.2 GeV) = 0.22). However, this term comes with a large 
logarithmic factor of ln(/z 2 /M^) = —6, for an appropriate scale of /1 = 4.2 GeV. The 
total correction to C± = 1 in ( fL9|) is then —21% (42%)! The presence of the large 
logarithm spoils the validity of a straightforward perturbative expansion, despite the fact 
that the coupling constant itself is still reasonably small. This situation is quite common in 
renormalizable quantum field theories. Logarithms appear naturally and can become very 
large when the problem involves very different scales. The general situation is indicated 
in the following table, where we display the form of the correction terms in higher orders, 




(19) 
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denoting I = \n(n/M w ) 



LL NLL 



OiJ 



2 1)2 



all 



a 



a s 

a\l 

all 2 



a: 



a: 



(20) 



a:: 



i i 
C(l) 0(a s ) 

In ordinary perturbation theory the expansion is organized according to powers of a s 
alone, corresponding to the rows in the above scheme. This approach is invalidated by 
the large logarithms since a s l, in contrast to a s , is no longer a small parameter, but 
a quantity of order 1. The problem can be resolved by resumming the terms (a s l) n 
to all orders n. The expansion is then reorganized in terms of columns of the above 
table. The first column is of 0(1) and yields the leading logarithmic approximation, 
the second column gives a correction of relative order a s , and so forth. Technically the 
reorganization is achieved by solving the renormalization group equation (RGE) for the 
Wilson coefficients. The RGE is a differential equation describing the change of C±(/j,) 
under a change of scale. To leading order this equation can be read off from ([19|) 

d 



d In [i ± ^ 4tt ^ ± 



C±(aO 



(21) 



(a s /47r)7±' ) are called the anomalous dimensions of C±. To understand the term "di- 
mension", compare with the following relation for the quantity /i" - , which has (energy) 
dimension n: 

d „ 



dlrifi 



n ■ \x 



(22) 



The analogy is obvious. Of course, the C±(fi) are dimensionless numbers in the usual 
sense; they can depend on the energy scale \x only because there is another scale, M w , 
present under the logarithm in (fl9|) . Their "dimension" is therefore more precisely called 
a scaling dimension, measuring the rate of change of C± with a changing scale /i. The 
nontrivial scaling dimension derives from 0(a s ) loop corrections and is thus a genuine 
quantum effect. Classically the coefficients are scale invariant, C± = 1. Whenever a 
symmetry that holds at the classical level is broken by quantum effects, we speak of an 
"anomaly". Hence, the 7±' ) represent the anomalous (scaling) dimensions of the Wilson 
coefficients. 



We can solve (|21|), using 

da <, 



and find 
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2/30 
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MAO ln M w 



(0) 
2/30 



(23) 



(24) 



4-7T ^ 

This is the solution for the Wilson coefficients C± in leading logarithmic approximation, 
that is to leading order in RG improved perturbation theory. The all-orders resummation 
of a s log terms is apparent in the final expression in (|24"D . 
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3.3 AB = 1 effective Hamiltonian 

In this section we will complete the discussion of the AB = 1 effective Hamiltonian. So 
far we have considered the operators 

Qi = {diPi)v-A{Pjbj)v-A (25) 
Ql = (dip^v-Alpjbijv-A (26) 

which come from the simple H^-exchange graph and the corresponding QCD corrections 
(Fig. §). We have slightly generalized our previous notation, allowing for the cases p = u, 
c. In addition, there is a further type of diagram at 0(a s ), which we have omitted until 




-u, c 
u, c 



Figure 8. QCD correction to W exchange. 
now: the QCD-penguin diagram shown in Fig. |9|. It gives rise to the four new operators 




Figure 9. QCD-penguin diagram. 



Q 3 = 


{dibi) V -A^{qjqj)v-A 
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(27) 


Q 4 = 


(dibj^AYsiWih-A 
q 


(28) 


Q 5 = 


{dibi) v ^ A Y^{mj)v+A 
q 


(29) 


Q 6 = 


(dib^v-Aj^iWiJv+A 


(30) 



q 



Two structures appear when the light-quark current (qq)v from the bottom end of the 
diagram is split into V — A and V + A parts. In turn, each of those comes in two colour 
forms in a way similar to Q\ and Q2- Finally, one further gauge-invariant operator of 
dimension six appears in the matching procedure, the chromomagnetic penguin operator 

Qs 9 = -^- 2 m b d t a^(l + l5 )T?^ G; v (31) 

This operator corresponds to the diagrams in Fig. |9| with the lower quark line omit- 
ted. The gluon is thus an external field, represented in ( |3~T1) by the field-strength tensor 
G a . Note that the characteristic tensor current necessitates a helicity flip in the b — ► d 
transition, which is accompanied by a factor of the quark mass m;, (the effect of is 
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neglected). The contribution of Qg g would be very small for the Hamiltonian of K decays, 
which only involves light external quarks, but it is unsuppressed for b decays. 

The operators Qi, ■ ■ ■ ,Qe, Qsg m i x under renormalization, that is the RGE for their 
Wilson coefficients is governed by a matrix of anomalous dimensions, generalizing (j2ip . 
In this way the RG evolution of Ci^ affects the evolution of C3, . . . , Cq, C$ g . On the other 
hand Ci j2 remain unchanged in the presence of the penguin operators Q 3 , . . . , Q e , Q 8g , so 
that the results for derived above are still valid. 

The construction of the effective Hamiltonian follows the principles we have discussed 
in the previous sections. First the Wilson coefficients Ci(fj,w), i = 1, • • -,6,8(7, are de- 
termined at a large scale fiw = 0(Mw,m t ) to a given order in perturbation theory. 
In this step both the W boson and the heavy top quark are integrated out. Since the 
renormalization scale is chosen to be nw = 0(M w ,m t ), no large logarithms appear and 
straightforward perturbation theory can be used for the matching calculation. The anoma- 
lous dimensions are computed from the divergent parts of the operator matrix elements, 
which correspond to the UV-renormalization of the Wilson coefficients. Solving the RGE 
the Ci are evolved from fiw to a scale /x = 0(rrib) in a theory with / = 5 active flavours 
q = u,d, s, c, b. The terms taken into account in the RG improved perturbative evaluation 
of Ci(n) are, schematically: 



LO: (a, In ^f) NLO: a s (a s In 



M w 



at leading and next-to-leading order, respectively. 

The final result for the AB = 1 effective Hamiltonian can be written as 



n eff ~ —ft 2^ A P 
V * p=u,c 



C\Q\ + C2Q2 + X] CiQi 

i=3,...,6,8s 



+ h.c. (32) 



where X p = V* d V p b. In principle there are three different CKM factors, X u , X c and X t , cor- 
responding to the different flavours of up-type quarks that can participate in the charged- 
current weak interaction. Using CKM unitarity, one of them can be eliminated. If we 



eliminate Aj, we arrive at the CKM structure of (|32|) . 

The Hamiltonian in (|32|) is the basis for computing nonleptonic 6-hadron decays within 
the standard model (to lowest order in electroweak interactions) with AB = 1 and AS, 
AC = 0. The Hamiltonian for 6-decays with different flavour quantum numbers of the 
light quarks has a completely analogous form. For instance AB = 1 transitions with 
a simultaneous change in strangeness, AS = 1, are simply described by (|32]) after the 
replacement d — > s. When new physics is present at some higher energy scale, the 
effective Hamiltonian can be derived in an analogous way. The matching calculation 
at the high scale nw will give new contributions to the coefficients Ci(nw), the initial 
conditions for the RG evolution. In general, new operators may also be induced. The 
Wilson coefficients Ci are known in the standard model at NLO. A more detailed account 
of 7~L^f B f =1 and information on the technical aspects of the necessary calculations can be 
found in (Buchalla, Buras & Lautenbacher 1996) and (Buras 1998). 
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3.4 B - B mixing at NLO 



In the following section we present the effective Hamiltonian for AB = 2 transition, which 
is relevant for B B mixing. In this case only a single operator contributes. The form 
of the Hamiltonian is therefore particularly simple. We use this example to illustrate the 
structure of Wilson coefficients at next-to- leading order. The mass difference Am^ in the 

/AB= 
L eff 



B - B system is related to the effective Hamiltonian H^M 2 by 



Am; 



2IM- 



12 | B 



1 

m B 



(B\H% B f 2 \B) 



(33) 



In order to construct Ti-^ff 2 ■• ^ ne standard model amplitude for AB = 2 transitions is 



matched onto the effective theory amplitude at the matching scale /i t 
This is sketched in Fig. [K]. 



0{m t ) = 0{M, 



w i 




b 

W +QCD ~ C(nt)- 

d 




d 



Figure 10. OPE for B - B mixing. 



There is only one local operator 

Q = {bd)v_ A {bd) v ^ A 
The Wilson coefficient, up to next-to- leading order, can be written as 



(34) 



(35) 



where is the lowest order result and C^ l > comes from the corrections with one-gluon 
exchange. The RG evolution from the high scale fj, t down to a scale /i = C?(ra&) has the 
form 

a s (n) - a s (ii tJ 

— J 5 



C(jm) 



1 + 



Arc 









_ « s (yU) . 



6/23 



(36) 



The second factor on the right-hand side is familiar from the leading logarithmic approx- 
imation (the only difference is that at NLO the two-loop expression for a s {fi) has to be 
used). The first factor represents the next-to-leading order correction. Here J5 is a scheme- 
dependent constant, which in the usual, so-called NDR scheme reads J 5 = 5165/3174. 

We now have the ingredients to write the effective Hamiltonian up to NLO precision 



n 



AB=2 
eff 



CQ*)Q 



r 2 m 2 

167T 2 ^ tb u) 



(37) 



1 . a 8 (fj) T 

1 + ^T Js 



Q 
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The result is entirely dominated by the top-quark contribution. It is common practice to 
separate the coefficient C(/x) into the function S (x t ) (x t = m 2 /M^), which would be the 
coefficient in the absence of QCD effects, into the terms that depend on a s (/i), and the 
remainder, which is defined as the QCD-correction factor rj B . This has been done in the 
second equation in fl37l) . Taking the matrix element of 7"^f/-/ =2 between the B and the B 
state and using (^) gives 



C 2 /If 2 

Am B = - J ^\V? b V td \ 2 S (x t ) VB B B f B m B 
One encounters the hadronic matrix element of Q, which is written as 



(3* 



(B\Q\B)fr) = -f*m 2 B B B (fi) 



(39) 



defining the (scale and scheme dependent) hadronic parameter B B {jj). The combination 

oc s (fi) 



B 



B 



B B (n)[a a (n))»a 



1 + 



Air 



-J, 



(40) 



is formally scale and scheme independent and has been used in fl38l) . The parameter B B {^i) 
is a nonperturbative quantity and has to be determined e.g. by lattice calculations. At 
present the value of B B is still very uncertain, in contrast to the short- distance QCD 
corrections, which are precisely known. A numerical illustration is given in Table 1, 
where we have put B B {^) = 0.9 as an example. 



Table 1. 



0.846 


0.9 


= 0.761 


rj B [a.(A*)]3 [l + ^^l BM 




0.551 


1.38 


= 0.761 



Two different definitions of a short- distance QCD factor can be considered, depending 
on where the terms with a s (fi) are included. One possibility is to include them with 
f] B into a Wilson coefficient (=0.846), which is to be multiplied by the hadronic matrix 
element B B {^) = 0.9. The other possibility is the formally scheme independent separation 
into r/ B = 0.551 and B B = 1.38 (for r\ B this is the precise result; B B = 1.38 is only true 
in our example). The purpose of this exercise is to remind us that different definitions 
are sometimes employed for the parameter B B and care has to be taken which one is 
being used, in order to combine it with the appropriate short-distance corrections. We 
can also see that the large deviation of the QCD correction factor r\ B from 1 is merely a 
consequence of pulling out the large factor [a s (/x)]~ 6//23 . It is somewhat artificial and does 
certainly not indicate a problem for perturbation theory. In fact, the coefficient 0.846 is 
the one that has the proper limit, approaching 1 as a s — > 0. It is indeed much closer to 
unity in accordance with the expectation for a perturbative correction factor. Still, the 
use of 7] B = 0.551 is often adopted due to its formally scheme invariant definition. 
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An important application is the ratio of the mass differences for and B s mesons, 
for which fl58"|) implies 



Am B d 



td 



V 



Is 



m B d f Bd B B d 

m B J 2 B B Bs 



(41) 



This quantity is a very useful measure of |V^/14s|. All other short- distance physics (top- 
dependence, i] B ) has dropped out. Hadronic uncertainties are reduced in the ratio of 
matrix elements, which is 1 in the limit of unbroken SU(3) flavour symmetry. The can- 
cellation of the short- distance contribution is a direct consequence of the factorization 
property of the OPE. Lattice calculations give for the ratio of matrix elements (Hocker 
et al. 2001, and refs. therein) 



1.16 ±0.06 



f Bd \ /B B 



(42) 



The ratio Am Bd / Am Bs is a very powerful constraint for the unitarity triangle, as can be 
seen in Fig. j|. 



4 Heavy quark effective theory 
4.1 Basic formalism 

Heavy quark effective theory is an effective field theory designed to systematically exploit 
the simplifications of QCD interactions in the heavy-quark limit for the case of hadrons 
containing a single heavy quark. The HQET Lagrangian can be derived as follows. We 
start from the usual QCD Lagrangian for a heavy-quark field \& with mass m 

£ = ^iJZ>^ - (43) 

with the covariant derivative 

D, = d, - igT a Al (44) 
The heavy-quark momentum can be decomposed as 

p = mv + k (45) 

where v is the 4- velocity of the heavy hadron. Once mv, the large kinematical part of the 
momentum is singled out, the remaining component k is determined by soft QCD bound 
state interactions, and thus k = 0(Aq C d) m. We next decompose the quark field \l/ 
into 

h v {x) = e mvx ^V(x) (46) 
H v (x) ee e mvx ^-^(x) (47) 

which implies 

^( x ) = e~ imv - x (h v (x) + H v {x)) (48) 

The expressions (1± )f)/2 are projection operators. Their action represents the covariant 
generalization of decomposing \1> into upper and lower components. Using the standard 
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representation for 7-matrices, this is evident in the rest frame where $ = 7 . Note also that 
the equation of motion with respect to the large momentum components, m(yt— l)h v = 0, 
is manifest for h„. 



The exponential factor exp(imv ■ x) in (|4q), fl47|) removes the large-frequency part of 
the x-dependence in $?(x) resulting from the large momentum mv. Consequently, the 
x-dependence of h v , H v is only governed by the small residual momentum and derivatives 
acting on h v and H v count as 0(Aq C d)- (Our sign conventions are appropriate for a 
heavy quark. To describe the case of a heavy anti-quark, similar definitions are valid with 
the sign of v reversed.) 

Multplying the QCD e.o.m. (i p — m)\I/ = with the projectors (1— fi)/2 and 
(1+ f)/2, and using (||) - @ and the definition 



D^ = D^- v»v ■ D (49) 

we obtain the coupled system of equations 

iv-Dh v = -%jp L E v (50) 
(iv ■ D + 2m)H v = ift L h v (51) 

They represent the e.o.m. in terms of h v and H v . The second equation implies that 
H v = 0(AQ C D/ m )h v by power counting. Hence H v is suppressed with respect to h v in 
the heavy-quark limit. In other words, h v contains the large components, H v the small 
components of 



The HQET Lagrangian is obtained starting from (f43|), expressing \1/ in terms of h v , 
H v and eliminating H v using d51|). We find 



C = h v iv ■ Dh v + h v i fl± i p±h v (52) 

iv ■ D + 2m 

Alternatively, H v as obtained from (|5T| ) in terms of h v can be inserted into (|50| ) to yield 
the e.o.m. for h v . This equation is just the e.o.m. implied by (^) (upon variation with 
respect to h v , i.e. 5C/5h v = 0). The Lagrangian may thus be written down immediately 
given the e.o.m. for the field h v . 

The second term in ( p2[ ) contains the nonlocal operator (iv ■ D + 2m)~ 1 . It can be 
expanded in powers of Aqcd/™ to yield a series of local operators. Keeping only the 
leading-power correction we can simply replace (iv ■ D + 2m)~ 1 by (2m) _1 and get 

C = h v iv ■ Dh v + —h v (iD ± ) 2 h v + —h v a^G, u h v (53) 
2m Am 

Let us discuss some important aspects of this result. 



• The first term on the r.h.s. of (^) is the basic, lowest-order Lagrangian of HQET. 
It describes the "residual" QCD dynamics of the heavy quark once the kinematic 
dependence on m is separated out. Since there is no longer any reference to the mass 
m, the only parameter to distinguish quark flavours, this term is flavour symmetric: 
The dynamics is the same for b and c quarks in the static limit. Since the operator 
v ■ D contains no 7-matrices, which would act on the spin degrees of freedom, the 
leading HQET Lagrangian also exhibits a spin symmetry. This corresponds to the 
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decoupling of the heavy-quark spin in the m — > oo limit. Together, we have the 
famous spin - flavour symmetries of HQET (Isgur & Wise 1989). They lead to 
relations among different heavy-hadron form factors. 

From the Lagrangian h v iv ■ Dh v the Feynman rules for HQET can be read off. The 
propagator is 

v ■ k 2 

where the projector (1+ $)/2 appears since h v is a constrained spinor (see fl46|) ). 
The interaction of the heavy-quark field h v with gluons is given by the vertex 

igvT" (55) 

These Feynman rules enter in the computation of QCD quantum corrections. 

The remaining terms in (|53|) are the leading power corrections. They have an 
intuitive interpretation. In the first term one recognizes the operator for the nonrel- 
ativistic kinetic energy / (2m) , which describes the residual motion of the heavy 
quark recoiling against the light degrees of freedom inside the heavy hadron. The 
last term represents the chromomagnetic interaction of the heavy-quark spin with 
the gluon cloud. Both effects violate flavour symmetry, the chromomagnetic term 
also spin symmetry, but they are power suppressed. 

So far we have only considered QCD interactions. Weak interactions introduce ex- 
ternal currents, which can also be incorporated in HQET. A generic heavy-light 
transition current qT^, arising for instance in semileptonic decays, can be repre- 
sented as 

qTV = qTh v + 0{— ) (56) 
m 

replacing the heavy-quark field \I/ by the HQET field h v using (|48|). 



4.2 Theory of heavy-hadron masses 

Before considering HQET in the context of weak decays, let us discuss a first application of 
the basic HQET Lagrangian ( p3[) in the spectroscopy of heavy hadrons. To be specific, we 
shall analyze the masses of the ground-state mesons B and B*. These mesons constitute 
a doublet that arises because the spin 1/2 of the heavy quark couples with the total spin 
1/2 of the light degrees of freedom in their ground state to form a spin-0 and a spin-1 
meson, the pseudoscalar B and the vector B*, respectively. Because the 6-quark spin 
decouples in the heavy-quark limit, the state of the light cloud is identical for B and B* 
to leading order, and the angular-momentum coupling described above is the appropriate 
scheme. If we neglect the power corrections in (0), we can immediately write down the 
composition of the meson masses 

to# = m^l = rrib + A (57) 

Evidently the meson mass has a component m& from the heavy quark. In addition it has 
a term A = O(Aqcd) from the energy of the light constituents. The latter is determined 
only by the interactions among the light degrees of freedom and their interaction with 
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the static 5-quark (h v ) through the first term in (p3|). It is therefore independent of m b . 
The sum of mj and A is a physical quantity, however, separately both parameters are 
dependent on the scheme used to define them. 

In order to include the first power corrections, we treat the 1/m terms in (|53| ) as 
perturbations to the lowest-order HQET dynamics. To first order in perturbation theory 



the corrections to (|57|) are then simply given by the expectation values of the 1/m terms. 
The proper normalization is obtained as follows. If TC — —Ci/ m is the Hamiltonian 
(density) corresponding to the correction term Ci/ m in ([53]), and H = J d 3 xTC is the 
Hamilton operator, the mass correction due to H is just 

5m B = (B^HlBi) (58) 

where \Bi) is the 5-meson state normalized to one, (Bi\Bi) = 1. Using the conventionally 
normalized states with (B\B) = 2m B V we can write 

5m B = ^y(B\ J d 3 xH(x)\B) = ^JL_ J d 3 x (B\H(0)\B) = {B ^ B) (59) 

where we have used the translation invariance of 7i and J d 3 x = V. Defining 

, _ (B\h{iDfh\B) 1 (B\hga ■ Gh\B) 

M = o A 2 = ^ z (60) 

IvriB o 2m B 

we obtain 

Sm B = — (61) 

Zmi, 

Note that we may replace D 2 by D\ in the definition of Ai, up to higher order corrections 
(see (^9]), (|50D). The parameter Ai corresponds to (minus) the expectation value of the 
momentum squared of the heavy quark, Ai = — (pf) = 0(Aq CD ). This gives a positive 
correction in Q5T|) representing the (small) kinetic energy of the heavy-quark. The A2- 
correction to the mass reflects the interaction energy of the heavy-quark spin with its 
hadronic environment, as already discussed in the previous section. While the Ai-term 
is independent of the heavy-quark spin and identical for B and B*, the chromo magnetic 
correction ~ A2 = 0(Aq CD ) is different for B*. We have 

Sm B * = (62) 

Including (|57|) we arrive at the following expansion for the meson masses 

1 t Ai + 3A 2 , s 

m B = m b + A (63) 

2m\ ) 

m B * = m b + A - — — (64) 

2m b 

where the dependence on m b is explicit order by order. 

If we apply the heavy-quark limit to D mesons, we obtain analogous relations 

■ t Ai + 3A 2 , , 

m D = m c + A — (65) 

m D * = m c + A- — — (66) 

2m r 
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with the same A, Ai and A2 as before. 

These results have a few interesting consequences. First, A2 parametrizes the spin- 
splitting between the pseudoscalar and the vector mesons: 

2A 2 

m B , - m B = —± = 46 MeV (67) 
m b 

2A 

m D * -m D = — - = 141 MeV (68) 
m c 

HQET predicts that the spin-splitting scales inversely proportional to the heavy-quark 
mass. This is seen to be quite well fulfilled given that m& ~ 3m c . Relation (|67D can be 
used to determine the nonperturbative quantity A2 from experiment 

A 2 = -(m|, - m 2 B ) = 0.12 GeV 2 (69) 

On the other hand, the quantity Ai has to be estimated theoretically. Finally one may 
introduce the spin-averaged masses 

_ m B + 3m B * j. Ai 

m B = = m b + A - - — (70) 

4 2mb 

m D + 3m D * j- Ai 

m D = i = m ^ k ~^ c W 

This eliminates A2 and yields the useful result 

m b -m c = (m B - m D ) \ 1 - l _ ) (72) 

\ 2m B m D J 

Since the Ai-correction is fairly small, the quark-mass difference is rather well determined, 
much better than individual quark masses. 



Exercise 

Derive the relative factor between the chrome-magnetic correction to the mass of the B 
and the B* meson. 

Solution: Denote the heavy-quark spin by s, the total spin of the light degrees of 
freedom by j and the total spin of the meson by J = s+j. The chromomagnetic 
field of the light cloud has to be proportional to j. Hence the energy of the interaction 
between this field and s is proportional to (s ■ j). Angular momentum algebra gives 
(2s -J) = J(J + 1) - s(s + 1) - j(j + 1), which is (-3/2) for B and (1/2) for B\ hence 
the relative factor (—1/3) of the A2-term in (|62|) with respect to (|8"TD. 



4.3 Heavy-light currents and fs 

The 5-meson decay constant f B is defined by the matrix element 

(0\A^\B{p)) = -if B m B v^ (73) 
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of the heavy-light axial vector current 

A„ = ^75* (74) 

Here q is the light-quark, \l/ the heavy-quark field in full QCD, with \P = b in the present 
case. The 5-meson momentum is p = wibv. 

Let us analyze in HQET, including QCD corrections. The expansion of in 
HQET to leading order in 1/m, but allowing for QCD effects, has the form 

A = C x {fi)A x + CMA 2 + O(-) (75) 

m 

M = qj^K A 2 = qufl&hv (76) 

The matching conditions at the 6-quark mass scale \x = m b are 

C 1 (m b ) = l + 0(a a ) C 2 (m b ) = 0(a s ) (77) 

To leading order in QCD only A x is present in HQET, with coefficient one. Radiative 
corrections at 0(a s ) modify C x and generate a new operator A 2 . Note that the matching 
calculation of the full-QCD current A onto HQET, leading to ([75|), is completely anal- 
ogous to the OPE procedure of constructing the effective weak Hamiltonian from the 
M^-exchange amplitude in the full standard model, which we have discussed in sec. 3. 
The difference is only that a 1/Myy expansion is performed in the latter case, and a l/m b 
expansion in the case of HQET. The basic philosophy is essentially the same. In partic- 
ular, a factorization of long and short-distance contributions is obtained: Contributions 
from large scales > /i, including the m^-dependence, is again contained in the coefficient 
functions C x>2 . Soft scales < fi are factorized into the hadronic matrix elements of A Xj2 . 

In contrast to the full-QCD current A, the HQET currents do have an anomalous 
dimension, reflecting a logarithmic dependence of fs on the heavy-quark mass at 0(a s ). 
The logarithms can be resummed by renormalization group methods, again in full analogy 
to the procedure in sec. 3. In leading logarithmic approximation (LLA) C 2 can be 
neglected and C x acquires the familiar form 



Ci(m) 



, \ -1 -2//9o 



(78 



Here the LLA assumes the hierarchy a s {m b ) 1, a s ln(m&//i) = 0(1), which holds in 
the heavy-quark limit {m b large, fi — 0(1 GeV)). 

To express fs in HQET via (|73|) , (|75|) and (^|), we need the matrix element of A x 

(0\A x \B(p)) = -i/(/i)Vrni^ (79) 

Since the dynamics of HQET is independent of mj, the reduced decay constant /(/i) is 
m b - independent. The only m 6 -dependence in (|79[) enters through a trivial factor v /me 
from the normalization of the 5-meson state, which in the usual convention is given by 

{B\B) = 2m B V (80) 
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Collecting the ingredients, ( ff5|) yields 



B 



v^b [ a s (fJ>) _ 



-2/Po 



(81) 



This expression for fs is true to leading order in the HQET expansion in KqcD/ m b and in 
leading logarithmic approximation in QCD. /(/x) in (|8~]J) is still a nonperturbative quantity 
to be determined by other methods. However, the dependence of fg on the heavy-quark 
mass is now explicit. Eq. (|8~ID implies the scaling behaviour fs ~ l/^ms, up to a 
calculable logarithmic dependence on m&. In principle such a relation can be used to 
relate to the analogous quantity /d for heavy mesons with charm. In practice, it turns 
out that the leading order scaling result for /# is not very well fulfilled even for the 6-mass 
scale and that subleading power corrections are important in this case. Nevertheless the 
result in (^l|) is of conceptual interest and can serve as a simple example of an application 
of HQET. 



4.4 Heavy-heavy currents: B — ► D^'Ip and 

One of the most important results of HQET is the extraction of V c b from exclusive semilep- 
tonic B — > D*W decay. We will here give a short outline of the main steps in this analysis. 
Starting point is the differential decay rate 

dT(B-,D'W) ^ 
aw 

in the kinematical variable w = v ■ v' where v and v' are the 4-velocities of B and 



D*, respectively. The dependence of fl82|) on \V c b\, the quantity of interest, is obvious, 
and K.(w) is a known kinematical function. Finally, T{w) contains the nontrivial QCD 
dynamics encoded in the B — > D* transition form factors. The corresponding matrix 
elements of the weak currents can be written in the heavy-quark limit as 

1 (D*{v',e)\&r^b\B{v)} = i^{w)e{^e,v\v) (83) 

1 (D*(v',e)\^ l5 b\B(v)) = eW f(l+%-(e-^;l (84) 
y /mD*'mB L J 

In the heavy-quark limit, that is to lowest order in HQET, all hadronic dynamics is 
expressed in a single function £(w), the Isgur-Wise function (Isgur & Wise 1989). In this 
limit we further have 

F(w) = ^(w) (85) 
Moreover, £ is absolutely normalized at the no-recoil point 

ai) = i (86) 

The no-recoil point w — 1 corresponds to the kinematical situation where the D* meson 
stays at rest in the rest frame of the decaying B (v' — v =^ w — 1). Measuring dT/dw at 
w = 1, \V c b\ can then be determined from (p2|) since all ingredients are known. Because 
w = 1 is at the edge of phase space, an extrapolation is necessary to find dT/dw\ w= i from 
the measured spectrum. 



Heavy Quark Theory 



23 



For a realistic analysis corrections to the heavy-quark limit need to be considered. An 
important property of B — > D*lv is that linear power corrections in HQET are absent, 
&i/m = 0, where m can be either m c or m b . Consequently the leading corrections enter 
only at second order and are thus greatly reduced. This result is known as Luke's theorem. 
The absence of linear corrections does not hold for B — > Dlv decays, hence the particular 
importance of B — > D*W. Including corrections, the lowest order approximation = 
£(1) = 1 is modified to 

•F(1) = tm(1 + «W) (87) 

where 5 1 / m 2 are the second order power corrections and tja is a correction from perturbative 
QCD. To first order in a s it reads 

VA = l + ^(^±^ln^-l) (88) 
7i Kirif, — m c m c 3/ 

A detailed numerical analysis yields (Harrison & Quinn 1998) 

^(1) =0.913 ±0.042 (89) 

which gives (Hocker et al. 2001) 

V cb = 0.0409 ± 0.00U exp ± 0.0019^ (90) 

To summarize the crucial points for the extraction of V cb from B — > D*W decay: 

• Heavy-quark symmetry relates the various semileptonic form factors (four different 
functions V, A Q , A 1 , A 2 in full QCD) to a single quantity £(u>), the Isgur-Wise 
function. 

• The function £ is absolutely normalized, £(1) = 1. This property has an intuitive 
reason: At the kinematical point w = 1 the decaying 6-quark at rest is transformed 
into a c-quark, also at rest. Since both quarks are treated in the static approximation 
(nib, trie oo, mb/m c fixed), the light hadronic cloud doesn't notice the flavour 
change b — > c and is transfered from the B to a D meson with probability one. The 
function £ is identical for B — > D and B — > D* transitions, because these are related 
by heavy-quark spin symmetry. 

• HQET provides a framework for systematic corrections to the strict heavy-quark 
limit governed by £(w). Luke's theorem guarantees the absence of first-order cor- 
rections in \ jm for B — > D*W. 



4.5 HQET - conclusions 

We would finally like to summarize the basic ideas and virtues of HQET, and to re- 
emphasize the salient points. 

• HQET describes the static approximation for a heavy quark, covariantly formu- 
lated as an effective field theory and allowing for a systematic inclusion of power 
corrections. 
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• Order by order in the expansion in Kqcd/ 171 HQET achieves a factorization of hard, 
perturbative contributions (momentum scales between m and a factorization scale /x) 
and soft, nonperturbative contributions (scales below //). The former are contained 
in Wilson coefficients, the latter in the matrix elements of HQET operators. 

• The procedure of matching full QCD onto HQET is analogous to the construction 
of the effective weak Hamiltonian H e ff- The difference lies in the massive degrees 
of freedom that are being integrated out: the W boson (mass M w ) for 7i. e ff, the 
lower-component spinor field H v (mass 2m) for HQET. The perturbative matching 
can be supplemented by RG resummation of logarithms, ln(Mw/fi) in the former 
case, ln(m///) in the latter. 

• The usefulness of HQET is based on two important features: The spin-flavour sym- 
metry of HQET relates form factors in the heavy-quark limit and thus reduces 
the number of unknown hadronic quantities. The dependence on the heavy-quark 
masses is made explicit (scaling, power corrections). 

We conclude with briefly mentioning another field, called large energy effective theory 
(LEET), which has some similarities with HQET. LEET is needed for matrix elements of 
the form (M\qT b\B) at large recoil of the light meson M = ir, p, K^*\ . . .. HQET is not 
sufficient in this situation because not only soft but also collinear infrared singularities 
need to be factorized. The latter occur due to the light-like kinematics of the fast and 
energetic light quark emitted from the weak current. To define LEET the usual heavy- 
quark limit can be considered for the B meson with velocity v. The large-energy limit is 
taken for the light meson M with light-like momentum vector En. Here E = 0{n r i} ) ) is 
the energy of M and n is a light-like 4- vector with n 2 = and v ■ n = 1. The momentum 
of the energetic light quark q is written as p q = En + k, with a residual momentum 
k = O(Aqcd)- In formal analogy to the fields h v and H v in HQET, the new light-quark 
fields 

q n (x) = e iE ™ A4(a;) Q n (x) = e iE ™ ^q(x) (91) 

can be defined and used in the construction of LEET (Dugan & Grinstein 1991, Charles 
et al. 1999, Beneke & Feldmann 2001, Bauer et al. 2001a). As a consequence of the 
LEET limit the ten form factors needed to describe all matrix elements (M\qT b\B) of 
bilinear heavy-light currents can be reduced to only three independent functions. LEET 
has received increasing interest quite recently and is still under active development. 

5 Inclusive decays and the heavy quark expansion 
5.1 Basic formalism and theory of lifetimes 

The heavy-quark limit, m ^> Aqcd, proves to be extremely useful also for the computation 
of inclusive decay rates of heavy hadrons (Chay et al. 1990, Bigi et al. 1992, 1997). The 
specific technique appropriate for this application is distinct from HQET and goes by 
the name of heavy quark expansion (HQE). Consider the total decay rate of a heavy 
hadron H. Starting point for the HQE is the following representation of T H 



(92) 
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where the transition operator T is defined as 

r.ta./AT^HH^ID) (93) 



with Tleff the effective weak Hamiltonian. Eqs. (p2|), (|93|) express the total decay rate 
as the absorptive part of the forward scattering amplitude H —>■ H under the action of 
Ti-eff- This expression is refered to as the optical theorem by analogy to a similar relation 
in optics. One may rewrite ©, (H) in a more directly understandable form by inserting 
a complete set of states between the two factors of 7~t e ff in ( F^D and removing the 

T-product by explicitly taking the absorptive part. This yields 

r H ~{H\?Uft\X){X\7Uff\H) (94) 

where one immediately recognizes the decay rate as the modulus squared of the decay 
amplitude (summed over all final states X). The reason to introduce fl9lf ) is that the T- 
product, by means of Wick's theorem, allows for a direct evaluation in terms of Feynman 
diagrams. 

In order to compute T H an operator product expansion is applied to (p3|), resulting in 
a series of local operators of increasing dimension. The coefficients of these operators are 
correspondingly suppressed by increasing powers of I /nib. The series has the form 

T = T b bb + % bga ■ Gb + V ^ bT t q qT i b + ... (95) 
mi m% 

where we have written the first few operators of dimension three (bb), five (bga ■ Gb) and 
six (bTiq qTib). The matrix elements of the operators contain the soft, nonperturbative 
physics, their Wilson coefficients r&, Zk the hard contributions, which are calculable in 
perturbation theory. Again, the coefficients are determined by an appropriate matching 
calculation between (1931) and the r.h.s. of (1951). The Feynman diagrams for the three 



terms in (R3) are shown in Fig. 11 The two weak-interaction vertices in these diagrams 
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Figure 11. Heavy quark expansion for the total decay rate of b-hadrons. 

correspond to the two factors of T~L e ff in the definition of T in ( p3|) (the absorptive part 
of the diagrams is understood). 
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Obviously, the heavy quark expansion is different from HQET. However, we may 
still use HQET in conjunction with (|95| ) in order to further analyse the hadronic matrix 
elements. An important example is the leading dimension-three operator bb. Its matrix 
element between heavy-hadron states H can be expanded in HQET as 

W = 1 + 7^71 CH*D) 2 h) + -r^ (hgaGh) (96) 



2ml Ami 



where (...) = (H\ . . . \H)/{2m H ). 



Eqs. (0), ( PBD and ( |9"E| ) imply that to leading order in the HQE T H = T b , that is the 
total decay rate of all 6-flavoured hadrons is equal to the rate of free 6-quark decay. 
Pictorially this can be seen from the first diagram in Fig. [LI], which represents 
essentially the amplitude squared for the partonic decay of a 6-quark. Note also 
that perturbative QCD corrections to r& can consistently be taken into account. 
The gluonic corrections to inclusive 6-quark decay are infrared safe, as required for 
Tfe in its role as a Wilson coefficient of the HQE. Also, corrections proportional to 
powers of a s (mi,) ~ 1/ m(m;,/A) are only suppressed by inverse powers of In m& in the 
heavy-quark limit, and hence formally leading in comparison to higher corrections 
in the HQE, which are suppressed by powers of Kjm b . The calculation of heavy- 
quark decay in the parton picture has been used since the beginnings of heavy-quark 
physics as an approximation for inclusive decays of the corresponding heavy hadrons. 
As we have seen, the HQE gives a formal justification for this approach and provides 
us with a theoretical framework to compute nonperturbative corrections. 



The first correction term in fl96|) depends on the expectation value of the momentum 
squared (p 2 ) of the heavy quark inside the hadron. This matrix element is non-zero 
because the heavy quark is recoiling against the light degrees of freedom through 
gluonic interactions in the hadronic bound state. This term has a very intuitive 
interpretation. It corresponds to a correction factor 1 — (p 2 )/(2m^) = 1 — (tf)/2, 
which is just the reduction of the free decay rate from time dilatation due to the 
recoil motion of the heavy quark. The second correction comes from interactions of 
the light hadronic cloud with the heavy-quark spin. We have 

f |(m 2 R , - ml) H = B 
(hgaGh) = B> „ . (97) 

{ H = A b 

The result is zero for the A& baryon since the light degrees of freedom are in a state 
of zero total angular momentum. Note that the spin interaction enters twice in ( |95*D , 
explicitly with coefficient zq and via the expansion of (bb). 

The leading nonperturbative corrections start only at second order. There is no 
correction linear in \jm b . This is because there is no gauge-invariant operator of 
dimension four that could appear in the HQE. 

At order \jm\ contributions appear where the spectator quark participates directly 
in the weak interactions. For 6-mesons they can be interpreted as the effect of weak 
annihilation of the 6-quark with the valence <i-quark (for Bo) and as the effect of 
Pauli interference (for B u ). The latter phenomenon occurs because in the nonlep- 
tonic decay of a B u , b{u) — > cud(u), two identical -u-quarks are present in the final 
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state. These corrections distinguish in particular among B d and B u mesons and are 
responsible for their lifetime difference. Despite the suppression by three powers 
of m b these effects can be relatively important due to their two-body kinematics, 
which brings a phase-space enhancement factor of 167T 2 in comparison to the leading 
three-body decay. 

As one of the possible applications, the HQE provides us with a theory of heavy- 
hadron lifetimes. The deviations of lifetime ratios from unity probes the power corrections. 
At present there are still sizeable theoretical uncertainties due to the hadronic matrix 
elements (bTqqTb). They can in principle be computed with the help of lattice gauge 
theory. Table 2 shows a comparison of theoretical predictions and experimental results 
(see for instance (Ligeti 2001)). 



Table 2. 





exp. 


th. 


t(B+)/t(B°) 


1.068 ±0.016 


1 - 1.1 


f(B s )/r(B d ) 


0.947 ±0.038 


0.99-1.01 


r{A b )/r{B d ) 


0.795 ±0.053 


0.9-1.0 



5.2 Local quark-hadron duality 

A systematic uncertainty within the HQE framework, which is often debated in the lit- 
erature, arises from the issue of quark-hadron duality. In this paragraph we give a brief 
and heuristic discussion of the basic idea behind this topic. 

The theoretical prediction for an inclusive decay rate obtained from the HQE has the 
form 

00 / A \ n 

T/T = l + Y.zn — (98) 

n=2 

where we have denoted the leading, free-quark result by r . Let us consider the decay 
rate as a function of m b , keeping A = Aq C d constant. Then the quantity r/r , to any 
finite order in (A/m), is a simple polynomial expression in this variable. This is sketched 
as the monotonous curve in Fig. |T2| showing r/r as function of m b (in arbitrary units). 
Now, since by construction the HQE for r/T yields a power expansion in (A/m), any 
term of the form 



exp I — I — I I sin 



^) 09) 



A J J V A ) 

for example, present in the true result for r/T would be missed by the HQE. This is 
due to the exponential suppression in the expansion parameter. In fact, the function 
exp(— 1/x) is non-analytic. Its power expansion around x = gives identically zero. 
However, such (or similar) terms are expected to be part of the true T/Tq on general 
grounds. The corresponding complete result for T/Tq, including such a term, is sketched 
as the oscillating graph in Fig. 0. This true curve represents the physical result for 
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Figure 12. r/r as function of m^. 

the decay rate T/Yq, which consists of the inclusive sum over all the different exclusive 
decay channels. It is intuitively understandable that the true m^-dependence will have 
such a damped oscillating behaviour: If we imagine to continually increase m?,, T /T$ will 
undergo a small jump whenever it reaches a value at which the presence of a further higher 
hadronic resonances in the final state becomes kinematically allowed. Since the excited 
hadrons have finite widths, the threshold behaviour will be smoothed out, resulting in the 
pattern of damped oscillations. 

The term quark-hadron duality refers to the idea that the inclusive rate as the sum 
over all exclusive hadronic decay channels and the inclusive rate as predicted by the heavy 
quark expansion are dual to each other. This means they are both valid representations 
of the same quantity using different descriptions, the hadron level or the quark level. The 
term local refers to the fact that the energy scale m& is a fixed quantity, as opposed to e.g. 
the centre-of-mass energy in e + e~ annihilation, which can be averaged to obtain suitably 
defined "global" quantities. In principle, the hadronic description gives the true result, 
measured in experiment. The problem is, however, that we would have to compute all 
exclusive rates first, which is far beyond our current control of nonperturbative QCD. 
On the other hand, the HQE calculation can be performed, within some uncertainties, 
but it is clear that the result need not be identical to the true answer. A deviation 
between the latter and the HQE (including power corrections) is refered to as a violation of 
quark-hadron duality. Indeed, contributions violating quark-hadron duality are expected 
(seefl99|)), but the numerical size of these terms cannot be strictly computed at present. 
Conceptually this is no problem because they are formally subleading in comparsion to 
power corrections, so that the HQE still makes sense even at higher orders. The remaining 
question is how large can violations of quark-hadron duality be numerically. While there 
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are at the moment, within the uncertainties intrinsic to HQE, no established cases in 
inclusive B decays where duality is violated, the issue clearly needs further investigation, 
both theoretically and phenomenologically. 

A more detailed account of the status of quark-hadron duality can be found in the 
papers by Blok, Shifman & Zhang (1998), Shifman (2000) and Bigi & Uraltsev (2001). 



5.3 Inclusive semileptonic decays — V u b, V c b 

The HQE cannot only be applied to the total decay rates, but also to inclusive rates 
with specific flavour quantum numbers in the final state, such as semileptonic processes. 
Furthermore one can analyze differential decay rates. 

An example of special interest is the inclusive decay B — ► X C W, which can be used to 
extract V cb . The HQE for the integrated rate has the form 



T(B^X c lu) = ^4\V cb \ 

K c 1 1927T 3 1 1 



Ai + 3A 2 \ 6A 2 



Z3 \ 1 + 2m? ) + Z5 ml + 



(100) 



with 



A -l- 3A - 3 

3 f 2ml 2 (btrCfb) = 6X2 = -(ml. - ml) (101) 



The Wilson coefficients read 

z 3 = 1 - 8x + 8a; 3 - a; 4 - 12x 2 In x + 0{a s ) (102) 
z 5 = -(1 - x) 4 (103) 

Where x = (m c /m(,) 2 . 

A major source of theoretical uncertainty for the determination of \V c b\ using (|100|) 
is the 6-quark mass. This appears to be especially problematic since m;, comes with the 
fifth power in (|100|) . Fortunately, however, the actual situation is not as bad. Taking into 
account the phase-space function z 3 , one finds that the combined dependence on m& and 
m c shows the approximate behaviour 

T(B -> X c lu) ~ m 2 b 3 (m b - m c f 7 (104) 

Since the difference m& — m c is better known than the individual quark masses, the 
corresponding uncertainty is reduced. The quark-mass difference is in fact constrained by 
HQET, which gives fl7|) 

m b -m c = (m B -m D )(l ^ 



2mBmo J 

= 3.40 ± 0.03 ± 0.03GeV (105) 

where tub = (tub + 3mB*)/4. 

The QCD corrections to are known to 0(a s ) and partly at 0(a 2 s ). The special class 
of corrctions (9(/3Q _1 a™) has been calculated to all orders n. 

Numerically the inclusive method gives (Hocker et al. 2001) 

V cb = 0.04076 ± 0.00050 e:C p ± 0.00204^ (106) 
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which can be compared with the result from the exclusive determination via B — > D*lv. 

One can also try to extract \V u b\ from B — > X u li> decays. This is more difficult since 
the very large background from semileptonic h — > c transitions requires kinematical cuts 
(in the lepton energy, the hadronic or the dilepton invariant mass), which renders the 
HQE less reliable and introduces larger uncertainties. A recent discussion has been given 
by Bauer et al. (2001b). The HQE has further useful applications, for instance in the 
case of the inclusive rare decays B — > X s ^y, B —>■ X s ^l + l~ , or B —>■ X sa -pp. 



Exercise 

Show that quark-hadron duality is exactly fulfilled for the semileptonic b — > c transition 
rate in the Shifman-Voloshin (small-velocity, or SV) limit rrib, m c 3> m& — m c ^> Aqcd- 
This holds with only two exclusive channels on the hadronic side of the duality relation, 
that is the inclusive rate is saturated as T(B — ► XJu) = T(B — > Dip) + T(B — > D*Ip) in 
this limit. 

Solution: We start from the exclusive differential decay rates in the heavy-quark limit. 
They read (see e.g. Harrison & Quinn 1998): 

dT(B^Dlp) _ Gl\V4\ {mB+mDfmUw2 _ lf/2e{w) {m 



dw 487r 3 
dT(B^D*lp) G 2 F \V cb \ 



■(m B — rn* D ) 2 rn z D ,\/w 2 — l(w + 1) 



2 



dw 487r 3 

( Aw m 2 B -2wm B m D ,+m 2 D A >} m 

In the strict SV limit we have 

m B = m h mp* = m D = m c m c = 7715(1 — e) (109) 

where e = (nib — n^c) / m b is a small parameter. 

The variable w is related to the dilepton invariant mass q 2 through 

q 2 = m 2 B + m 2 D — 2m B m D w (HO) 

The kinematic limits of q 2 are easily identified as 

(Lax = ( m b ~ m c ) 2 q 2 min = (111) 

The corresponding limits of w are 

m 2 + m 2 

yJmin 1 ^ 'max — ~ (H^) 

lmbm c 

Defining s = w — 1 we have < s < e 2 /2, where the upper limit is valid to leading order 
in e. Expanded to leading order in e, ( |107| ) gives 

r* 2 n/ i 2 /-e 2 /2 r^ 2 it/ 1 2 

r(S - DZv) = ^fUnf^ / , 3 / 2 rf, = ^Mm 6 - m c ) 5 (113) 

07T' 5 JO O07T J 
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which is the decay rate in the SV limit. In this derivation we have made use of the fact 
that £(u>) = £(1) + C(e 2 ), which can be approximated by £(1) = 1. In this way any 
dependence on nontrivial hadronic input has disappeared. Similarly we can expand the 
integral over ( |108| ) in e to extract the leading contribution in the SV limit. We obtain 

F(B -> D*lv) = ^f^K - ™ c ) 5 (H4) 

We also observe that higher D-meson resonances and hadronic multiparticle states have 
wave functions of the light degrees of freedom that are orthogonal to the ground state 
wave function of the light cloud (identical for D, D* and B) in the SV limit. There is 
therefore no overlap of those higher excitations with the initial B and the corresponding 
rates vanish. 

Finally, we need to take the SV limit of the inclusive rate as obtained from the heavy 
quark expansion in ( |100| ). In this limit the second-order power corrections and perturba- 
tive QCD corrections disappear, and we only have to expand the phase space function z 3 
in the small-e limit. We find z% = 64e 5 /5 + 0(e 6 ) and 

T(B - XJu) = G \^ c f (m b - m c f (115) 
loir-* 

We see that indeed the inclusive HQE result (|H5|) is saturated by the sum of just the 
two exclusive rates (|113|) and (|114j) . Clearly, the SV limit is a very special situation. 



Nevertheless, it is an interesting example of exact (local) quark-hadron duality. Moreover, 
the semileptonic rates into D and D* measured in experiment account for roughly two 
thirds of the inclusive rate, indicating that the SV limit is not even entirely unrealistic. 



6 QCD factorization in exclusive hadronic B decays 
6.1 Introduction 

Decay amplitudes for exclusive nonleptonic B decays, such as B — > ttti, can be computed 
starting from the effective weak Hamiltonian discussed in sec. 3.3. Whereas the Wilson 
coefficients C, are well understood, the main problem is posed by the hadronic matrix 
elements of the operators Qi. In some cases this problem can be circumvented (CP 
asymmetry in B —>■ J/^fKs), or at least reduced using SU(2) or SU(3) flavour symmetries 
and an appropriate combination of various channels. However, an improved understanding 
of the QCD dynamics in exclusive hadronic B decays would greatly enhance our capability 
to extract from these processes the underlying flavour physics. 

Indeed, it turns out that the heavy-quark limit leads to substantial simplifications also 
in the problem of hadronic two-body decays of heavy hadrons. Again the main feature 
is the factorization of short- distance and long-distance contributions. In the case of the 
matrix elements of four-quark operators Qi the factorization takes the form 

(7r(p>(g)|Q,|B(p)) = f B ^(q 2 ) f duTl{u)^{u) 

Jo 

+ / dtdudvTf'fau^BiO&M&^v) (116) 
Jo 
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This factorization formula is valid up to corrections of relative order KqcD/ m b- Here 
f B ^ n (q 2 ) is a B — > n form factor evaluated at q 2 = m 2 ps 0, and $„. (<&b) are leading- 
twist light-cone distribution amplitudes "wave functions") of the pion (B meson). These 
objects contain the long-distance dynamics. The short-distance physics, dominated by 
scales of order m.&, is described by the hard-scattering kernels T?' 11 , which are calculable 
in perturbation theory. T/ starts at 0(a®), T/ 7 at 0{a\) (see Fig. |l3|). In (|116|) long- 




Figure 13. Order a s corrections to the hard scattering kernels T? (first two rows) and 
T/ 7 (last row). In the case of Tf , the spectator quark does not participate in the hard 
interaction and is not drawn. The two lines directed upwards represent the two quarks 
forming the emitted pion. 

and short- distance contributions are thus systematically disentangled, that is factorized. 
The long-distance sensitive quantities (form factors and wave functions) still need to be 
determined by other means, but they are universal quantities and much simpler than the 
original full B — > tttt matrix elements we started with. They could in principle be cal- 
culated by nonperturbative methods or extracted experimentally from other observables. 
In any case (|116 ) represents a substantial simplification of our problem. 



The general expression (|116|) further simplifies when we neglect perturbative a s - 
corrections. The T 11 term is then absent and the kernel T 1 becomes a constant in u, 
such that the pion distribution amplitude integrates to the pion decay constant. The 
matrix element of operator Q±, for instance, reduces to 

(Tr + n-\{vb) V - A (du)v-A\B) - (Tr + \{ub) v . A \B) • (7r|(du)y_ A |0) = im 2 B f B ^{Q)U (117) 

This procedure, termed "naive factorization" has long been used in phenomenological 
application, but the justification had been less clear. An obvious issue is the scheme 
and scale dependence of the matrix elements of four-quark operators, which is needed to 
cancel the corresponding dependence in the Wilson coefficients. This dependence is lost 
in naive factorization as the factorized currents are scheme independent objects. In QCD 
factorization ( |116|) the proper scale and scheme dependence is recovered by the inclusion 
of 0(a s ) corrections as we will see explicitly below. 
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A qualitative justification for (|117|) had been given by Bjorken (1989). It is based on 
the colour transparency of the hadronic environment for the highly energetic pion emitted 
in B decay (the tt~ in the above example, which is being created from the vacuum). This 
is related to the decoupling of soft gluons from the small-size colour- singlet object that 
the emitted pion represents. The QCD factorization approach as encoded in ( |116| ) may be 
viewed as a consistent formalization and generalization of Bjorken's colour transparency 
argument. This treatment of hadronic B decays is based on the analysis of Feynman dia- 
grams in the heavy quark limit, utilizing consistent power counting to identify the leading 
contributions. The framework is very similar in spirit to more conventional applications 
of perturbative QCD in exclusive hadronic processes with a large momentum transfer, as 
the pion electromagnetic form factor (see the article by Sterman and Stoler (1997) for a 
recent review). It justifies and extends the ansatz of naive factorization. In particular 
the method includes, for B — > tttt, the hard nonfactorizable spectator interactions, pen- 
guin contributions and rescattering effects (Fig. [13]). As a corollary, one finds that strong 
rescattering phases are either of 0(a s ), and calculable, or power suppressed. In any case 
they vanish therefore in the heavy quark limit. QCD factorization is valid for cases where 
the emitted particle (the meson created from the vacuum in the weak process, as opposed 
to the one that absorbs the 6-quark spectator) is a small size colour- singlet object, e.g. 
either a fast light meson (tt, g, K, K*) or a J/^f. 

Note that the term factorization is used here for two a priori entirely different things: 
In the case of QCD factorization ( |lltj| ), it refers to the factorization of short- distance 
and long-distance contributions. In the sense of the phenomenological approach of naive 
factorization (|1 1 7f) , it simply denotes the separation of the hadronic matrix element of 
a four-quark operator into two factors of matrix elements of bilinear currents. It is a 
nontrivial result that the latter, naive factorization is obtained as the lowest order ap- 
proximation of QCD factorization. To avoid confusion, it is useful to keep the distinction 
in mind. For example, the hard gluon exchange corrections between the two quark cur- 
rents in Fig. [13] are "nonfactorizable" in the sense of naive factorization, although they 
are a consistent ingredient of ( |116| ), hence "factorizable" in the sense of QCD. 

In the following we shall discuss QCD factorization in some detail using the example 
of B — ► Dtt decays. In this case the b — > u transition current is replaced by a heavy-heavy 
b — > c current. This case is somewhat simpler than B — > tttt since the spectator interaction 
(the T 11 term, bottom line of Fig. [13]) does not contribute to leading power. This is 
because for a heavy-to-heavy transition the spectator quark, and hence the gluon attached 
to it, is always soft. This leads to a suppression, according to the colour transparency 
argument, when this gluon couples to the emitted pion. Also penguin contributions are 
absent for B — > Dir. We shall illustrate explicitly how factorization emerges at the one- 
loop order in this specific case, and in the heavy-quark limit, defined as m&, m c 3> Aqcd 
with m c /rrib fixed. 

Further details on QCD factorization in B decays and additional literature can be 
found in the articles by Beneke et al. (1999, 2000 & 2001). 
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6.2 B —> Dtt: Factorization to one-loop order 
6.2.1 Preliminaries 

The effective Hamiltonian relevant for B — ► Dti can be written as 

Gi 



with the operators 



O = cry" (I -7 5 )6rf7 M (l -7b)m, 

8 = 07^(1-75)^6^(1-75)^. 



;il8) 



(119) 
(120) 



Here we have chosen to write the two independent operators in the singlet-octet basis, 
which is most convenient for our purposes, rather than in the more conventional bases 
of Qi, Q 2 or Q + , Q_ (see the discussion in sec. 3; because all four quark flavours are 
different in (|118|) , penguin operators are absent). The Wilson coefficients C , C$ have 
been calculated at next-to-leading order in renormalization-group improved perturbation 
theory (Altarelli et al. 1981, Buras & Weisz 1990) and are given by 



Co 



Nc + l 

2N r 



C+ + 



2N r 



where 



~a s (M w y 


d± - 







1 + 



4tt 

a s (M w 



a s (n) 



47T 



{B ± - J±) 



(121) 

(122) 
(123) 



(The coefficients C , C$ are related to the ones of the standard basis by C = G\ + C 2 /3 
and C 8 = 2C 2 .) We employ the next-to-leading order expression for the running coupling, 



An 



/? ln(/x 2 /A^cD) 



UN -2f 



Pi = 3 {N? - f Nc f - 2C F f 



(3, \n\n(^/A 2 QCD ) 
(51 ln(/. 2 /A 2 QCD ) 



2N r ' 



where N c is the number of colours, and / the number of light flavours. Aqcd 
the QCD scale in the MS scheme with / flavours. Next we have 



(o) 
7± 

2/V 



(o) 
7± 



±12 



2N r ' 



B 



± 



' 2N r 



B. 



(124) 

(125) 
A (/) is 

(126) 



The general definition of J± may be found in (Buchalla, Buras & Lautenbacher 1996). 
Numerically, for N c = 3 and / = 5 



_6_ 
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1587' 



(127) 
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(a) (b) (c) 

Figure 14. Basic quark-level topologies for B — > £>7r decays (q = u, d): (a) class-I, (b) class- 
II, (c) weak annihilation. Bd — > D + ir~ receives contributions from (a) and (c), Bd — > D°tt° 
from (b) and (c), and B~ — > D°tt~ from (a) and (b). Only (a) contributes in the heavy-quark 
limit. 



The quantities @o, Pi, d±, B± — J± are scheme independent. The scheme dependence of the 
coefficients at next-to- leading order is parametrized by B± in ( 122 ). In the naive dimen- 
sional regularization (NDR) and 't Hooft-Veltman (HV) schemes, this scheme dependence 
is expressed in a single number B with -Bndr = 11 and -Bhv — 7. The dependence of 
the Wilson coefficients on the renormalization scheme and scale is cancelled by a corre- 
sponding scale and scheme dependence of the hadronic matrix elements of the operators 
Oq and 0$. 

Before continuing with a discussion of these matrix elements, it is useful to consider 
the flavour structure for the various contributions to B — > Dn decays. The possible quark- 
level topologies are depicted in Fig. |14[ In the terminology generally adopted for two-body 
non-leptonic decays, the decays Bd — > D + n~, Bd — > D n and B~ — > D°tt~ are referred 
to as class-I, class-II and class-Ill, respectively. In both Bd — > D + n~ and B~ — > D°7r~ 
decays the pion can be directly created from the weak current. We may call this a class-I 
contribution, following the above terminology. In addition, in the case of Bd —>■ D + n~ 
there is a contribution from weak annihilation and a class-II amplitude contributes to 
B~ — > D°7T~, see Fig. [14|. The important point is that the spectator quark goes into 
the light meson in the case of the class-II amplitude. This amplitude is suppressed in 
the heavy-quark limit, as is the annihilation amplitude. The amplitude for Bd — > D°7i°, 
receiving only class-II and annihilation contributions, is therefore subleading compared 
with Bd — > D + n~ and B~ — > D°tt~, which are dominated by the class-I topology. The 
treatment of this leading class-I mechanism will be the main subject of the following 
sections. We shall use the one-loop analysis for Bd — > D + tt~ as a concrete example on 
which we will illustrate explicitly how the factorization formula can be derived. 



6.2.2 Soft and collinear cancellations at one-loop 

In order to demonstrate the property of factorization for Bd — > D + tt~, we have to an- 
alyze the "non-factorizable" one-gluon exchange contributions (Fig. TE) to the b — > cud 



transition. We consider the leading, valence Fock state of the emitted pion. This is justi- 
fied since higher Fock components only give power-suppressed contributions to the decay 
amplitude in the heavy-quark limit. The valence Fock state of the pion can be written as 

Kg)) = / ^^VW C ( a T&) 6 ito) - a %) b \( 1 *)) |0> (128) 
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(a) (b) (c) (d) 

Figure 15. "Non-factorizable" vertex corrections. 

where a\ (b\) denotes the creation operator for a quark (antiquark) in a state with spin 
s =| or s =1, and we have suppressed colour indices. This representation of the pion 
state is adequate for a leading-power analysis. The wave function *&(u, l±) is defined 
as the amplitude for the pion to be composed of two on-shell quarks, characterized by 
longitudinal momentum fraction u and transverse momentum The on-shell momenta 
(lq,q = 0) °f the quark (l q ) and the antiquark (l q ) are given by 

I 2 I 2 
l q = uq + l L + —n_, l q = uq-l ± + ^n_. (129) 

Here q = E(l, 0, 0, 1) is the pion momentum, E = ps ■ q/rriB the pion energy and n_ = 
(1,0,0,-1). Furthermore l± ■ q = l± ■ n_ = 0. For the purpose of power counting 
lj_ ~ Aqcd <C E ~ m&. Note that the invariant mass of the valence state is (l q + lq) 2 = 
l]_/(uu), which is of order Aq CD and hence negligible in the heavy-quark limit, unless u 
is in the vicinity of the endpoints (0 or 1). In this case the invariant mass of the quark- 
antiquark pair becomes large and the valence Fock state is no longer a valid representation 
of the pion. However, in the heavy-quark limit the dominant contributions to the decay 
amplitude come from configurations where both partons are hard (u and u both of order 
1) and the two-particle Fock state yields a consistent description. The suppression of the 
soft regions (u or u 1) is related to the endpoint behaviour of the pion wave function. 
We will provide an explicit consistency check of this important feature later on. 

As a next step we write down the amplitude 

(n(q)\u(0)J(y)p\0) = J du^-^=^*(u,l ± )( l5 ^ a pe il ^, (130) 

which appears as an ingredient of the B — > Dtt matrix element. The right-hand side of 
(|130|) follows directly from (|128|) . Using (|130|) it is straightforward to write down the one- 



gluon exchange contribution to the B — ► Dtt matrix element of the operator 0% (Fig. |15|) 
We have 

(D + TT-IOglAiH-gluon = (I3i; 

2Cf f d*k t .. A . 1 ri 



where 



Am ~ 2 Vc -k-k 2 2 Pb -k + k 2 ' {U2) 
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AM I M_ r(^+#)7A lx(K+W 
A 2 {l q ,l q ,k)- 2l _. k + k2 2lq . k + k 2- 

Here T = 7^(1—75) and pi,, p c are the momenta of the b quark and the c quark, respectively. 
Note that this expression holds in an arbitrary covariant gauge. The gauge-parameter 
dependent part of the gluon propagator gives no contribution to (|131|) , as can be seen 
from (11321) and (11331). There is no correction to the matrix element of Oq at order a s , 



because in this case the (du) pair is necessarily in a colour-octet configuration and cannot 
form a pion. 

In ( |131| ) the pion wave function \l/(u, l±) appears separated from the B — > D transition. 
This is merely a reflection of the fact that we have represented the pion state by ( |128| ). 
It does not, by itself, imply factorization, since the right-hand side of ( |131| ) involves 
still nontrivial integrations over l± and gluon momentum k, and long- and short- distance 
contributions are not yet disentangled. In order for ( |131|) to make sense we need to 
show that the integral over k receives only subdominant contributions from the region of 
small k 2 . This is equivalent to showing that the integral over k does not contain infrared 
divergences at leading power in l/m^. 

To demonstrate infrared finiteness of the one-loop integral 

J = Jd 4 k^ Mk) ® A 2 (l q , l q , k) (134) 

at leading power, the heavy-quark limit and the corresponding large light-cone momentum 
of the pion are again essential. First note that when k is of order mb, J ~ 1 for dimensional 
reasons. Potential infrared divergences could arise when k is soft or when k is collinear 
to the pion momentum q. We need to show that the contributions from these regions are 
power suppressed in m?,. (Note that we do not need to show that J is infrared finite. It 
is enough that logarithmic divergences have coefficients that are power suppressed.) 

We treat the soft region of integration first. Here all components of k become small 
simultaneously, which we describe by scaling k ~ A. Counting powers of A (d 4 k ~ A 4 , 
1/k 2 ~ A~ 2 , 1/p ■ k ~ A -1 ) reveals that each of the four diagrams (corresponding to 
the four terms in the product in (|134|) ) is logarithmically divergent. However, because 
k is small, the integrand can be simplified. For instance, the second term in A 2 can be 
approximated as 



7 a(^+ m Jx(u 4+h + ^ i-+ j£)T q x 



2l q -k + k 2 2uq ■ k + 2l± ■ k + ■ k + k 2 1' k 



(135) 



where we used that $ to the extreme left or right of an expression gives zero due to the 
on-shell condition for the external quark lines. We get exactly the same expression but 
with an opposite sign from the other term in A 2 and hence the soft divergence cancels 
out when adding the two terms in A 2 . More precisely, we find that the integral is infrared 
finite in the soft region when l± is neglected. When l± is not neglected, there is a diver- 
gence from soft k which is proportional to l^/m 2 . ~ Aq CD /7?t|. In either case, the soft 
contribution to J is of order AQ CD /m b or smaller and hence suppressed relative to the 
hard contribution. This corresponds to the standard soft cancellation mechanism, which 
is a technical manifestation of colour transparency. 
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Each of the four terms in ( |134j ) is also divergent when k becomes collinear with the 
light-cone momentum q. It is convenient to introduce, for any four-vector v, the light-cone 
components 



v± 



J°±V 3 



V2 



The collinear region is then described by the scaling 



k: 



A°, k\ ~ A, k' 



(136) 



(137) 



Then d 4 k ~ dk + dk~d 2 k ± ~ A 4 and q ■ k = q + k~ ~ A 2 , k 2 = 2k+k~ + k\ ~ A 2 . The 
divergence is again logarithmic and it is thus sufficient to consider the leading behaviour 
in the collinear limit. Writing k = aq + . . . we can now simplify the second term of as 



2/(7 ' /l ~T~ /u^ 



2(w + a)r 



(138) 



No simplification occurs in the denominator (in particular l± cannot be neglected), but the 
important point is that the leading-power contribution is proportional to q\. Therefore, 
substituting k = aq into A\ and using q 2 = 0, we obtain 



q\Ai 



2ap c ■ q 



2ap b ■ q 



0. 



:i39) 



employing the equations of motion for the heavy quarks. Hence the collinearly divergent 
region is seen to cancel out via the standard collinear Ward identity. This completes 
the proof of the absence of infrared divergences at leading power in the hard-scattering 
kernel for Bd — > D + n~ to one-loop order. In other words, we have shown that the "non- 
factorizable" diagrams of Fig. ^ are dominated by hard gluon exchange. The proof at 
two loops has been given by Beneke et al. (2000) and a proof to all orders by Bauer et 
al. (2001c). 

Since we have now established that the leading contribution to J arises from k of order 
nib ("hard" k), and since E, we may expand A 2 in \l±\/E. To leading power the 

expansion simply reduces to neglecting l± altogether, which implies l q = uq and lg = uq 
in (|133|) . As a consequence, we may perform the l±_ integration in ( 131|) over the pion 



wave function. Defining 



d 2 l ± V*(u,lj 



167T 3 y/2W c 

the matrix element of Os in ( |131| ) becomes 



if* 



{D+K-\O s \B d ) X - 



gluon 



(140) 



:i4r 



'^W I 4^^ D+ ^ Al ^ Bd) h fn Jo du ®M tr i^dMuqiUq,k)}. 



Putting y on the light-cone in ( p.30| ), y + = y± = 0, hence l q ■ y = lty~ = uqy, we see 
that the /^-integrated wave function ^(u) in (|140|) is precisely the light-cone distribu- 
tion amplitude of the pion in the usual definition. Indeed, the leading-twist light-cone 
distribution amplitude for pseudoscalar mesons (P) reads 



(P(q)\q(y) a q'(xh\0) 



ifi 



(x-y) 2 =0 



5) 0a 



[ l du e i{Uqx+uqy) $j 
Jo 



(142) 



Heavy Quark Theory 



39 



Here it is understood that the operator on the left-hand side is a colour singlet. We 
use the "bar" -notation, i.e. v = 1 — v for any longitudinal momentum fraction variable. 
The parameter \x is the renormalization scale of the light-cone operators on the left-hand 
side. The distribution amplitude is normalized as Jq du $p(w, /i) = 1. One defines the 
asymptotic distribution amplitude as the limit in which the renormalization scale is sent 
to infinity. The asymptotic form is 

$ P (u,fi) M =°°6mx. (143) 

The decay constant appearing in ( |142| ) refers to the normalization in which f„ = 131 MeV. 
There is a path-ordered exponential that connects the two quark fields at different posi- 
tions and makes the light-cone operators gauge invariant. In (|142|) we have suppressed 
this standard factor. 

This derivation demonstrates the relevance of the light-cone wave function for the 
factorization formula. Note that the collinear approximation for the quark- ant iquark 
momenta emerges automatically in the heavy-quark limit. 

After the k integral is performed, the expression ( |141| ) can be cast into the form 



(£> + 7r-|0 8 | J B d ) 1 _ gluon ~ F B ^ D (0) f 1 duT 8 (u,z)^(u), (144) 

where z = m c /mb, T 8 (u,z) is the hard-scattering kernel, and Fb->d(Q) the form factor 
that parametrizes the (_D + |c[. . ]b\Bd) matrix element. The result for T 8 (u, z) will be given 
in the following section. 



6.2.3 Matrix elements at next-to-leading order 

As we have seen above, the — > D + n~ amplitude factorizes in the heavy-quark limit 
into a matrix element of the form (D + |c[. . ]b\Bd) for the B — > D transition and a matrix 
element (7r~|d(x)[. . .]w(0)|0) with x 2 = that gives rise to the pion light-cone distribution 
amplitude. Leaving aside power-suppressed contributions, the essential requirement for 
this conclusion was the absence of both soft and collinear infrared divergences in the gluon 
exchange between the (cb) and (du) currents. This gluon exchange is therefore calculable 
in QCD perturbation theory. We now present these corrections explicitly to order a s . 

The effective Hamiltonian ( |118| ) can be written as 



G p 

V2 



2N C vr/ 4vr 2N c 



+ C 8 (fi)0 8 \, (145) 



where the scheme-dependent terms in the coefficient of the operator On, proportional to 
the constant B defined after (|127|) , have been written explicitly. 

Schematically, the matrix elements of both On and 8 can be expressed in a form 
analogous to (|144j ). Because the light-quark pair has to be in a colour singlet to produce 
the pion in the leading Fock state, only On gives a contribution to zeroth order in a s . 
Similarly, to first order in a s only 8 can contribute. The result of computing the diagrams 
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in Fig. 15 with an insertion of Og can be written in a form that holds simultaneously for 
H = D,D* and L = ir,p, using only that the (ud) pair is a colour singlet and that the 
external quarks can be taken on-shell. One obtains (z = m c /m^) 



(H(p')L(q)\0 8 \B d (p)) 



g s C F 
An 2N r 



if L / du$ L (u) 



(146) 



x 



6 In 



mi 



+ B ){{J V )- (J a )) + F(u, z) (J v ) - F(u, -z) {J A ) 



where 



(J v ) = (H{p')\c4b\B d {jp)), (J A ) = (H(p')\c^ 5 b\B d (p)). (147) 

In obtaining (|146|) we have used the equations of motion for the quarks to reduce the 
operator basis to Jy and J a- The form of ( |146| ) is identical for pions and longitudi- 
nally polarized p mesons. The production of transversely polarized p mesons is power 
suppressed in Aq CD /m fo . 

In the case of a distribution amplitude $l(m) that is symmetric under u <-> u, which 
is relevant for L = tc, p, the function F(u, z) appearing in ( 146 ) can be compactly written 

as 

(148) 

with 



f(u, z) 



F{u,z) = 3\nz 2 - 7 + f(u,z) + f(u, l/z), 
u(l - z 2 )[3(l - u(l - z 2 )) + z] 



[1 - u{\ - z 2 )} 2 



\n[u(l-z 2 )} 



l-u(l- z 2 )' 



(149) 



The contribution of f(u, z) in (|148| ) comes from the first two diagrams in Fig. [15] with the 
gluon coupling to the b quark, whereas f(u, l/z) arises from the last two diagrams with 
the gluon coupling to the charm quark. The absorptive part of the amplitude, which is 
responsible for the occurrence of strong rescattering phases, arises from f(u, l/z) and can 
be obtained by recalling that z 2 is z 2 — ie with e > infinitesimal. We then have 



7T 



lm F(u, z) 



(1 - u)(l - z 2 )[3(l - u(l - z 2 )) + z] 
[l-u(l-z 2 )} 2 



(150) 



As mentioned above, ( |146| ) is applicable to all decays of the type B d — > D^ + L^, 
where L is a light hadron such as a pion or a (longitudinally polarized) p meson. Only 
the operator Jy contributes to B d — > D + L~ , and only Ja contributes to B d — > D* + L~ . 
(Due to helicity conservation the vector current B — > D* matrix element contributes only 
in conjunction with a transversely polarized p meson and hence is power suppressed in 
the heavy-quark limit.) Our final result can therefore be written as 



(D+L-\0 Q , s \B d ) = (D+\crf(l-<y 5 )b\B d )-if Lqil [* duT , 8 (u, z) <Z> L (u) 

Jo 

where L = it, p, and the hard-scattering kernels are 
T (u,z) 



Ta(ll, Z) 



l + 0(a 2 s ), 



6 In 



if! 

ml 



B + F(u,z) 



+ 0(al 



(151) 

(152) 
(153) 
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When the D meson is replaced by a D* meson, the result is identical except that F(u, z) 
in Q153Q must be replaced by F(u, —z). Since no order ot s corrections exist for O , the 
matrix element retains its leading-order factorized form 



(D+L-\O \B d ) = if LQfl {D + \vf(l - js)b\B d ) 



(154) 



to this accuracy. From ( |149|) it follows that T$(u, z) tends to a constant as u approaches 
the endpoints iu — > 0, 1). Therefore the contribution to (|151|) from the endpoint region 
is suppressed, both by phase space and by the endpoint suppression intrinsic to $l(m), 
which can be represented as 



<f> L (u) = 6u(l 



u 



n=l 



(155) 



vanishing ~ u (~ (1 — u)) at the endpoints. Here C^ 2 (x) = 3x, C*2 /2 (x) = |(5a; 2 — 1), 
etc. are Gegenbauer polynomials. The Gegenbauer moments a^(/i) are nonpertrubative 
quantities. They are multiplicatively renormalized and approach zero as fi — > oo. (The 
scale dependence of these quantities enters the results for the coefficients only at order 



,3/2, 



or 



which is beyond the accuracy of a NLO calculation.) 



As a consequence of the endpoint suppression the emitted light meson is indeed dom- 
inated by energetic constituents, as required for the self-consistency of the factorization 
formula ( |151| ). 

Combining ( |145| ), ( |151[ ), ( |152| ) and ( |15ij| ), we obtain our final result for the class-I, non- 
leptonic — ■> Lr^> + L~ decay amplitudes in the heavy-quark limit, and at next-to- leading 
order in a s . The results can be compactly expressed in terms of the matrix elements of a 
"transition operator" 

Gp 



T 



V2 



V* d V cb [ai(DL) Q v - a x {D*L) Q A \ 



(156) 



where 



Qv = cy% <g> d'j^l - 7 5 )m, Q a = 07^756 ® ^(1 - 7 5 )u, 



(157) 

and hadronic matrix elements of Qv,a are understood to be evaluated in factorized form, 
i.e. 

(DL\j x ®j 2 \B) = (D\j x \B) (L|j 2 |0). (158) 

Eq. (|156|) defines the quantities a x (D^L), which include the leading "non-factorizable" 
corrections, in a renormalization-scale and -scheme independent way. To leading power in 
these quantities should not be interpreted as phenomeno logical parameters (as 
is usually done), because they are dominated by hard gluon exchange and thus calculable 
in QCD. At next-to- leading order we get 



a x (DL) 



+ 



-61n^ + 



duF(u, z) 



(159) 



ax(D*L) 



N c + 1 
2N r 



CM + 



Nr. 



2N r 



, as Cp n ( \ 



u 2 r 1 
-61nAr+/ duF(u, -z)$ L (u) 
mi Jo 



(160) 
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These expressions generalize the well-known leading-order formula 



We observe that the scheme dependence, parametrized by B, is cancelled between the co- 
efficient of Oq i n Q145 ) and the matrix element of Og in ( |151| ). Likewise, the /i dependence 



of the terms in brackets in ( |159| ) and ( 160 ) cancels against the scale dependence of the 
coefficients C±(fi), ensuring a consistent physical result at next-to-leading order in QCD. 

The coefficients a\{DL) and a\(D*L) are seen to be non-universal, i.e. they are explic- 
itly dependent on the nature of the final-state mesons. This dependence enters via the 
light-cone distribution amplitude of the light emission meson and via the analytic 

form of the hard-scattering kernel (F(u,z) vs. F(u, —z)). However, the non-universality 
enters only at next-to-leading order. 



Exercise 

Verify that the coefficients ai(D^L) in ( \15£\ ) and ( \1 6(\) are independent of the unphysical 
renormalization scale /i up to terms of 0(a 2 s ). 



6.2.4 Phenomenological applications for B — > Dtt 

An important test of QCD factorization can be performed by comparing the hadronic 
decays B d — > D^ + L^ with the semileptonic processes B d -> DM+rz/. It is useful to 
introduce the ratios 



where X p = X* = 1 for a vector meson (because the production of the lepton pair via 
a V — A current in semi-leptonic decays is kinematically equivalent to that of a vector 
meson with momentum q), whereas X v and X* deviate from 1 only by (calculable) terms 
of order m 2 /m 2 B , which numerically are below 1%. The main virtue of ( |162| ) is that 
the B — > form factors cancel in the ratio. The theoretical prediction for the QCD 
coefficients, based on QCD factorization to leading power and at NLO in perturbative 
QCD, is \ai(D^L)\ = 1.05. The uncertainty of this leading-power result is small, about 
±0.01. The prediction is to be compared with the experimental results, extracted from 



( EH) , which read |ai(D*7r)| = 1.08 ± 0.07, |ai(D*p)l = 1.09 ± 0.10 and |a x (D*ai)| = 
1.08 ± 0.11. These values show fair agreement with the theoretical number, albeit within 
experimental uncertainties that are still large. 

Another interesting consideration concerns the comparison of class-I modes with those 
of class II and III. For B — > D^tt, all three decay modes, which are related by isospin, have 
been measured. A nice discussion of the present experimental status and its interpretation 
in the context of QCD factorization has been given by Neubert & Petrov (2001). Let us 
briefly discuss here a few important aspects. The experimental status is summarized 
in Table 3. Denoting the basic topologies from Figs. (a), (b) and (c), by T, C 



Heavy Quark Theory 



43 



Table 3. Experimental data for B — ► D^n branching ratios (in units of 10 3 ) and 
extracted parameters \C — A|/|T + A\, 5tc (see Neubert & Petrov (2001)). 







5 -> D*7T 


B- _> D(*)°tt- 


3.0 ±0.4 
0.27 ±0.05 
5.3 ±0.5 


2.76 ±0.21 
0.17 ±0.05 
4.6 ±0.4 


|C-A|/|T + A| 


0.42 ±0.05 
(56 ±20)° 


0.35 ±0.05 
(51 ±20)° 



and A, respectively (where the notation refers to "tree", "colour-suppressed tree" and 
"annihilation"), we have 

A(B d -> D + n-) = T + A (163) 

V2A{B d -> flV) = C-4 (164) 

^(fi-^D°7r-) = T + C (165) 

For later use we may further define the spectator-interaction contribution to T, T spec (see 
bottom row of Fig. [13]). A similar decomposition holds for the B — > D*tt modes. Isospin 



°7r°) 



symmetry is reflected in the amplitude relation A(B d — > D + n~) + ^/2A(B d — > D 
A(B~ — > D°7T~), which is manifest in the parametrization of (|163|) - (|165|) . This means 
that there are only two independent amplitudes, which we can take to be (T + A) and 
(C — A). These amplitudes are complex due to strong phases from final-state interactions. 
Only the relative phase of the two independent amplitudes is an observable and we define 
5tc to be the relative phase of (T + A) and (C — A). The phase can then be extracted 
from the data through the relation 

4§Xb(B~ -> D°ti-) - B(B° -> D+Ti-) - 2B(B° -> D tt°) 
cos S T c = — — ; _ _ (166) 



2^^JB(B° -»• D + 7T-) 5(5° -»• L>°7r°) 



In the usual heavy-quark limit, m b ~ m c ^> Aq C d, only T is calculable. T spec , C and 
A are not, but they are power suppressed. It is instructive to consider the alternative 
limit mj ^> m c 3> Aqcd, which allows us to distinguish between mj, and m c (Beneke et 
al. 2000). In this case, due to m c <C mj,, also the D becomes a "light" meson and in 
that respect the process is similar to B — > tttt, where both T spec and C are calculable. 
However, since m c ^> Aqcd = A, the D-meson wave function $d(m) is highly asymmetric 
and strongly peaks at u = (1 — w) ~ A/m c , where u is the c-quark momentum fraction. 
These properties can be used to derive the scaling rules 



A ^ A_ Tspec ^A_ 

T mi, T m c T m c 



(167) 



The amplitude A is still not calculable in this scheme, while T, T spec and C are. Note 
that from ( |167| ) we can recover the two standard scenarios we have been discussing: In 
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the heavy-quark limit m c ~ m& ( |167| ) reduces to a simple power suppression ~ A /m& for 
A, T spec and C compared to T. On the other hand, if m c becomes a truly light quark, 
corresponding to the case of B — * tttt, we count m c ~ A and see that both T spec and C 
are of the same order as T, while A is still power suppressed. 

The general scenario m& ^> m c ^> Aqcd allows us to interpret the experimental 
results in Table 3. We can even make numerical estimates for T, T spec and C based on 
the factorization formula for light-light final states (|116|) . These are somewhat model 
dependent because $£> is not known at present. It is not too difficult to accomodate 
substantial values |C — A\/\T + A\ ~ 0.2 - 0.3 and Stc ~ 40°, in qualitative agreement 
with Table 3. Given the special role of the charm quark (not light, but also not too 
heavy), the current experimental situation is not in contradiction with QCD factorization 
in the large-m^ limit. For a comparison with experiment it is useful to keep in mind that, 
according to (|167|) , the suppression of C over T is only ~ A/m c (not A/m b ) and that Stc 



can also be substantial. 
Exercise 

Derive the relation ( 166] ). 



6.3 CP violation in B — ► tt + tt decay 

Hadronic B decays into a pair of light mesons, such as B — ► ttK or B — > tttt, have a very 
rich phenomenology. Their main interest lies in their sensitivity to short-distance flavour 
physics, including CKM parameters, CP violation and the search for new physics. By 
way of an outlook we mention here the important example of CP violation in B — > tt+tt- 
decay. The starting point for computing the required decay amplitudes is the effective 
Hamiltonian in fl32|) . The needed matrix elements of the operators Qi can be analyzed 
within QCD factorization using ( |116| ). We will not go into the technical details of such 
an analysis and the discussion of limitations of the approach, in particular from power 
corrections in AgcD/n^b- These can be found in (Beneke et al. 2001). Here we just want 
to present the phenomenological motivation and to illustrate that a theoretical approach 
towards a direct dynamical calculation of hadronic matrix elements will be very valuable, 
even if it necessitates some approximations. 

The observable of interest is the time-dependent CP asymmetry in the decays B°,B° — > 
7r + 7r~, which is sensitive to the B&-B& mixing phase e~ 2tf3 . We define 



A^(t) 



B{B°{t) -> TT + TT') - B{B°{t) -> TT + TT') 
B(B°(t) -> TT + TT-) + B(B°(t) -> TT + TT') 

—Skk sin(Am B t) + C W7T cos(Arne t) (168) 



where 



2 Im X 7T7T 1 | Avi-tt | 2 _ -2i/3 e n + Ptttt/T^ (IRQ} 

+ p' + 12' A ^~ e e^ + P IT 1 ' 

* K v 7T7r I - 1 - * I /X 7T7T I ° 1 ± 7T7T / 7T7T 

The amplitudes T n7T ("tree") and P nn ("penguin") are the components of the B — > tt+tt- 
amplitude corresponding to the terms in (^) involving \ u and A c , respectively. In the 
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standard phase conventions A c is real and X u has a weak phase —7, which has been factored 
out above. The coefficient C n7T , which is a function of 7, represents direct CP violation 
and is expected to be small. We shall not discuss it further here. The mixing-induced 
asymmetry SV-rr depends on 7 and /3. In fact, in the limit where P^/T^ is set to zero it 
follows that A^tt = e _2 ^ /3+7 - ) = e 2ia , and hence S n7T = sin 2a. In this limit X n7T is just the 
relative weak phase between the direct amplitude B — > 7t + tt~ and the one with mixing 
B — > B — ► 7r + 7r~. All dependence on hadronic input has canceled in this situation. In 
practice, however, P 7rn /T n7T is not fully negligible. It is here that information on hadronic 
dynamics becomes crucial. QCD factorization predicts that P n7T /T n7T is suppressed (either 
by a s or by powers of AQCD/mb), because T n7T can arise at tree level, P nn only through 
loops. Estimates within this framework give values of about 0.25 - 0.3. 
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Figure 16. Relation between sin 2a and the mixing-induced CP asymmetry S wn , assuming 
sin 2/3 = 0.48. The dark band reflects parameter variations of the first kind, the light band 
shows the total theoretical uncertainty. The lower portion of the band refers to values 
45° < a < 135°, the upper one to < a < 45° (right branch) or 135° < a < 180° — (3 
(left branch). 



To illustrate the effect of penguins, we first assume that |Kift/Vd,| and the weak phase 
(3 have been determined accurately. Then using 7 = 180° — a — (3 the expression for 
A,,-,,- in ( |169| ) becomes a function of a and our prediction for the penguin-to-tree ratio 
P 7Tn /T 7Tn . If we further assume that the unitarity triangle lies in the upper half of the 
(p, fj) plane, then a measurement of S n7T determines sin 2a with at most a two-fold discrete 



ambiguity. Figure |16| shows the relation between the two quantities for the particular case 
where IV^/V^I = 0.085 and (3 = 14.3°, corresponding to sin 2(3 = 0.48. The dark band 
shows the theoretical uncertainty due to input parameter variations, whereas the light 
band indicates the total theoretical uncertainty including estimates of the effect of power 
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corrections. We observe that for negative values sin 2a as preferred by the global analysis 
of the unitarity triangle, a measurement of the coefficient S n7T could be used to determine 
sin 2a with a theoretical uncertainty of about ±0.1. Interestingly, for such values of sin 2a 
the "penguin pollution" effect enhances the value of the mixing-induced CP asymmetry, 
yielding values of S W7r between —0.5 and —1. Such a large asymmetry should be relatively 
easy to observe experimentally. 




Figure 17. Allowed regions in the (p,fj) plane corresponding to constant values of the 
mixing-induced asymmetry S nn assuming the Standard Model. The widths of the bands 
reflect the total theoretical uncertainty. The corresponding bands for positive values of 
S n7T are obtained by a reflection about the p axis. The light circled area in the left-hand 
plot shows the allowed region obtained from the standard global fit of the unitarity triangle 
(Hocker et al. 2001). 



Although it illustrates nicely the effect of "penguin pollution" on the determination 
of sin 2a, Figure is not the most appropriate way to display the constraint on the 
unitarity triangle implied by a measurement of S n7T . In general, there is a four- fold discrete 
ambiguity in the determination of sin 2a, which we have reduced to a two-fold ambiguity 
by assuming that the triangle lies in the upper half-plane. Next, and more importantly, we 
have assumed that iV^/Vy and (3 are known with precision, whereas a is undetermined. 
However, in the Standard Model and the angles a, (3, 7 are all functions of the 

Wolfenstein parameters p and f\. It is thus more appropriate to represent the constraint 
implied by a measurement of S nn as a band in the (p, ff) plane. To this end, we write 

T i 7 = PT if} _ 2l0 = (1 - P) 2 -ff- 2ifj(l - p) P n7r = r„ e**™ 
' "v^+? (l-p) 2 + rf VF+rF 1 ' 

where r n7T e^ 1 * is independent of p and fj. We now insert these relations into (|169|) and 
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draw contours of constant S nn in the (p, fj) plane. The result is shown by the bands in 
Figure |T7[ The widths of the bands reflect the total theoretical uncertainty (including 
power corrections). For clarity we show only bands for negative values of S n7T ; those 
corresponding to positive S nn values can be obtained by a reflection about the p axis (i.e., 
fj — > — fj) . Note that even a rough measurement of S n7T would translate into a rather narrow 
band in the (p,fj) plane, which intersects the ring representing the \V u b/V c b\ constraint at 
almost right angle. In a similar way, the constraint is also quite robust against hadronic 
uncertainties. Even the approximate knowledge of hadronic matrix elements, as provided 
by QCD factorization, will therefore be very valuable and can lead to powerful constraints 
on the Wolfenstein parameters. 

7 Summary 

In these lectures we have discussed the theory of heavy quarks, focussing on the important 
case of B physics. (All methods relying on the heavy-quark limit could in principle be 
applied to charmed hadrons as well, but they are in most cases much less reliable due to 
the smaller value of the charm mass.) We shall conclude by summarizing the key points. 

• A crucial and general idea for dealing with QCD effects is the factorization of short- 
distance and long-distance dynamics. We have encountered this principle in many 
different forms and applications: 

— The OPE to construct the effective weak Hamiltonians (W^ -1 ' 2 ) factorizes 
the short- distance scales of order Mw, m t from the scales of order mj. 

— The heavy-quark scale m treated as a short- distance scale can be factorized 
further from the intrinsic long-distance scale of QCD, Aqqd- This leads to 
a systematic expansion of observables simultaneously in 1/m and a s (m) with 
often very important simplifications. The precise formulation of this class of 
factorization depends on the physical situation and can take the form of HQET, 
LEET, HQE or QCD factorization in exclusive hadronic B decays. 

• HQET exhibits the spin-flavour symmetry of QCD in the heavy-quark limit, which 
allows us to relate different form factors, and makes the tuq dependence explicit. 
Examples of typical applications are B -> D^lv or f B . 

• HQE is a theory for inclusive B decays. It justifies the "parton model" and allows 
us to study the nonperturbative power corrections. This is of great use for processes 
as B — > X c ^ u lv, B — > X s 7, B — > X s l + l~ , and the lifetimes of 6-flavoured hadrons. 

• QCD factorization, finally, refers to a framework for analyzing exclusive hadronic 
B decays with a fast light meson as for instance B — > Dn, B — > tttt, B — > ttK and 
B -> V-f. 

With these tools at hand we are in a good position to make full use of the rich experi- 
mental results in the physics of heavy flavours. We can determine fundamental parameters 
of the flavour sector, such as V u b, V c b, V t d/V ts , r\ and sin 2a, and probe electroweak dy- 
namics at the quantum level through b — > s^y or B - B mixing. This will enable us to 
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thoroughly test the standard model and to learn about new structures and phenomena 
that are yet to be discovered. 
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